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Preface 


The  goal  of  this  research  was  to  establish  a  means  of  designing  and  measuring  Ver¬ 
tical  Cavity  Surface  Emitting  Lasers  (VCSELs).  The  design  described  and  methods  of 
characterization  used  herein  constitute  the  first  capability  to  build  and  measure  these  de¬ 
vices  at  the  Air  Force  Institute  of  Technology,  and  Wright  Patterson  AFB.  The  capabilities 
for  fabrication,  and  tools  for  more  sophisticated  measurement  practices,  will  build  on  the 
lessons  learned  from  this  effort.  I  want  to  express  sincere  thanks  to  several  people  and 
organizations  for  their  help  in  this  research.  I’d  like  to  thank  my  thesis  advisor,  Capt 
Jeff  Grantham,  for  his  willingness  to  share  his  expertise  in  this  area;  especially  the  lab 
configurations,  where  I  had  very  little  experience.  His  first  hand  knowledge  of  VCSEL 
experimentation  was  irreplaceable  and  allowed  the  accomplishment  of  twice  my  expecta¬ 
tions.  Mr  Rick  Patton  deserves  my  appreciation  also.  His  technical  expertise  and  ability 
to  get  what  the  experiment  needed  on  a  moment’s  notice  was  remarkable.  I’d  also  like  to 
thank  Mr.  Greg  Smith  for  his  help  with  the  dye  laser,  the  attenuator,  and  the  Ti-Sapphire 
system,  which  saved  me  days  at  least  (weeks, . .  .maybe). 

Thanks  also  to  Wright  Labs/ERL,  who  grew  the  back  mirror  perfectly,  RADC/- 
OCPA,  who  sponsored  the  work  and  provided  funds  for  several  pieces  of  critical  equipment, 
and  the  University  of  Arizona,  who  provided  several  VCSEL  wafer  samples  and  growth 
information. 

Finally,  I  thank  my  wife  Yvonne  and  children,  J.C.,  Marcus,  and  Kristina,  who 
provided  the  love  and  support  that  was  necessary  to  carry  me  through  this  effort.  Although 
I  know  the  kids  were  good,  Yvonne  still  had  her  hands  full.  I  owe  her  tremendously.  You 
can  start  your  vacation  tomorrow,  Yvonne. 


Richard  J.  Bagnell 


Table  of  Contents 


Page 


Preface .  ii 

Table  of  Contents  .  iii 

List  of  Figures .  vii 

Abstract .  x 

I.  Introduction .  1-1 

1.1  Background .  1-1 

1.1.1  DBR  Mirrors .  1-2 

1.1.2  VCSEL  Cavity .  1-3 

1.1.3  VCSEL  Characteristics .  1-4 

1.2  Problem  Statement .  1-6 

1.3  Scope .  1-6 

1.4  Approach/ Methodology .  1-7 

1.5  Sequence  of  Presentation  .  1-7 

II.  Theory .  2-1 

2.1  Efficient  Vertical  Cavity  Design .  2-1 

2.2  Cavity  Design .  2-3 

2.2.1  F:  Confinement  Factor .  2-3 

2.2.2  Cavity  Losses .  2-5 

2.2.3  Design  Equations .  2-6 

2.3  Active  Region  Design  .  2-9 

2.3.1  Bulk  Active  Regions .  2-9 

2.3.2  Quantum  Well  Active  Regions .  2-12 


iii 


Page 

2.4  Mirror  Design .  2-19 

2.4.1  Transmission  Matrix .  2-19 

2.4.2  Refractive  Indices  .  2-23 

2.4.3  Reflectivity  Calculations .  2-24 

2.4.4  Spectral  Reflectivity  Stop  Band .  2-25 

2.4.5  Penetration  Depth .  2-25 

III.  Application  and  Performance .  3-1 

3.1  Introduction  .  3-1 

3.2  MQW  InGaAs  950  nm  Vertical  Cavity  Laser .  3-1 

3.2.1  Active  Medium .  3-1 

3.2.2  InGaAs  Active  Region  Parameters .  3-2 

3.2.3  WeU  Width  and  x  Concentration .  3-3 

3.2.4  Cavity  Construction .  3-5 

3.2.5  Mirror  Design .  3-6 

3.3  Carrier  Guided  VCSELs .  3-11 

3.4  Injection  Lasers:  Electrical  Pumping .  3-12 

3.4.1  Active  Region  Design  .  3-13 

3.4.2  DBR  Mirror  Design .  3-13 

rV.  Experimental  Method .  4-1 

4.1  Laboratory  Configuration .  4-1 

4.1.1  Dye  Laser  Configuration  Reflectivity  Measurements  4-5 

4.1.2  Additional  Reflectivity  Measurements .  4-5 

V.  Experimental  Results .  5-1 

5.1  Measurement  Parameters .  5-1 

5.1.1  Reflectivity .  5-1 

5.1.2  Lasing  Wavelength .  5-2 

iv 


Page 

5.1.3  Pump  Wavelength .  5-2 

5.1.4  Polarization .  5-2 

5.1.5  Mode .  5-2 

5.2  Sample  1:  3  A  UGS  GaAs  Active  Region .  5-3 

5.2.1  Material  Parameters .  5-3 

5.2.2  Sample  1  Measurements .  5-4 

5.2.3  Sample  1  Analysis .  5-8 

5.3  Sample  2:  3  A  UGS  GaAs  Active  Region .  5-10 

5.3.1  Material  Parameters .  5-10 

5.3.2  Sample  2  Measurements .  5-11 

5.3.3  Sample  2  Analysis .  5-15 

5.4  Design  VCSEL  Back  Mirror .  5-18 

5.4.1  Material  Parameters .  5-18 

5.4.2  Back  Mirror  Measurements .  5-18 

VI.  Conclusion  .  6-1 

6.1  Characterization .  6-1 

6.1.1  Spectral  Reflectivity .  6-1 

6.1.2  Polarization .  6-3 

6.1.3  Laboratory  Configuration  Improvements .  6-4 

6.2  Design .  6-4 

Appendix  A.  Spectral  Reflectivity  FORTRAN  Programs .  A-1 

A.l  Program  to  Calculate  the  Spectral  Reflectivity  in  Fabry-Perot 
Etalon  with  Distributed  Bragg  Reflector  Mirrors  (VCSEL  Struc¬ 
ture)  .  A-1 

A.2  Distributed  Bragg  Reflector  Mirror  Spectral  Reflectivity  Calcu¬ 
lation  Program .  A-8 

Bibliography .  BIB-1 


V 


List  of  Figures 


Figure  Page 

1.1.  Vertical  Cavity  Surface  Emitting  Laser  (VCSEL  Geometry)  .  1-2 

1.2.  Periodic  Gain  Structure  (PGS)  Cavity  .  1-4 

1.3.  Fabry  Perot  Effect  in  Typical  VCSEL .  1-5 

2.1.  Periodic  Gain  Structure  (PGS)  VCSEL  Cavity  Geometry .  2-4 

2.2.  Gain  vs  Nominal  Current  Density  Curve  for  GaAs  .  2-10 

2.3.  g/J„om  vs  Curvp  for  GaAs .  2-11 

2.4.  JthlJoiR)  vs  Fill  Factor(d//)  Curve  for  GaAs  .  2-13 

2.5.  Quantum  Well  Energy  Diagram .  2-14 

2.6.  Compressive  Strain  Effects  on  Valence  Sub-Bands  (GaAs) .  2-18 

2.7.  Geometry  for  Reflected  Fields  in  Multilayer  Stack .  2-21 

2.8.  Reflection  Delay  in  Multilayer  Stack .  2-26 

3.1.  Required  Quaintum  Well  Width  vs  Indium  Concentration,  A  =  950  nm  .  3-4 

3.2.  Cavity  for  the  Ino.2Gao,8As/GaAs  64  AQuantum  Well  VCSEL .  3-6 

3.3.  Bottom  Mirror  Reflectivity  as  seen  in  Air  .  3-7 

3.4.  Bottom  Mirror  Reflectivity  as  seen  in  GaAs  Spacer  .  3-8 

3.5.  Top  Mirror  (22  layers)  Reflectivity  as  seen  from  GaAs  Cavity  Spacer  .  3-10 

3.6.  Fabry- Perot  Effect  in  Designed  950  nm  VCSEL .  3-11 

3.7.  Gain  Guiding  and  Index  Guiding  in  VCSEL  Structures .  3-11 

3.8.  Gain  and  Index  Guiding  from  Carrier  Injection  .  3-12 

3.9.  Methods  to  Avoid  High  Mirror  Resistances  for  Current  Injected  VCSELs  3-14 

4.1.  Equipment  Configuration  used  for  Output  Beam  Characterization  ...  4-1 

4.2.  Optical  Pump  Beam  and  Focusing  Lens  Geometry  .  4-2 

4.3.  Additional  Reflectivity  Measurement  Configuration .  4-6 

vii 


Figure  Page 

5.1.  Sample  1:  3  A  Active  Region  VCSEL  Geometry .  5-3 

5.2.  Sample  1:  Theoretical  Spectral  Reflectivity  near  875  nm  .  5-4 

5.3.  Sample  1  Measured  Spectral  Reflectivities  for  Lasing  and  Non-Lasing  Lo¬ 
cations  near  Spot  1  .  5-5 

5.4.  Sample  1  Measured  Spectral  Reflectivities  for  Lasing  and  Non-Lasing  Lo¬ 
cations  near  Spot  2  .  5-6 

5.5.  Sample  1  Measured  Spectral  Reflectivity  for  Non-Lasing  Locations:  Spots 

5,  6,  and  7 .  5-7 

5.6.  Sample  1,  Spot  1  :  TEMqo  Lasing  and  Pump  Spots .  5-8 

5.7.  Sample  2:  3  A  Active  Region  VCSEL  Geometry .  5-10 

5.8.  Sample  2:  Theoretical  Spectral  Reflectivity  near  875  nm  .  5-11 

5.9.  Sample  2,  Spot  3:  TEMio  Lasing  and  Pump  Spots  .  5-12 

5.10.  Sample  2  Measured  Spectral  Reflectivities  for  Lasing  and  Non-Lasing  Lo¬ 
cations  near  Spot  3  .  5-12 

5.11.  Sample  2,  Spot  4:  TEMqo  Lasing  and  Pump  Spots  .  5-13 

5.12.  Sample  2  Measured  Spectral  Reflectivities  for  Lasing  and  Non-Lasing  Lo¬ 
cations  near  Spot  4  .  5-14 

5.13.  Sample  2  Measured  Spectral  Reflectivity  for  Non-L«ising  Locations:  Spots 

5  and  6 .  5-15 

5.14.  Sample  2  Measured  vs  Theoretical  (875  nm  VCSEL)  Spectral  Reflectivi¬ 
ties  for  Lasing  and  Non-Lasing  Areas .  5-16 

5.15.  Sample  2  Measured  Spectral  Reflectivities  for  Lasing  and  Non-Lasing  Ar¬ 
eas  vs  Theoretical  (887  nm  GaAs  VCSEL) .  5-16 

5.16.  Spectral  Reflectivity  Measurements  at  Center  of  Designed  920  nm  Back 

Mirror  .  5-19 

5.17.  Spectral  Reflectivity  Measurements  at  10  mm  radial  out  from  Center  of 

Designed  920  nm  Back  Mirror .  5-20 

5.18.  Spectral  Reflectivity  Measurements  at  Center,  10mm  radial,  and  20  mm 

radial  of  Designed  920  nm  Back  Mirror .  5-21 


viii 


Page 


Figure 

5.19.  Spectral  Reflectivity  Measurements  at  two  radials,  5mm  from  Center  of 


Designed  920  nm  Back  Mirror .  5-22 

6.1.  Non-Coincident  Pump  and  Exit  Beam  Locations  due  to  Surface  Rough¬ 
ness  .  6-3 


IX 


AFIT/GE/ENP/92D-01 


Abstract 

Vertical  Cavity  Surface  Emitting  Lasers  (VCSELs)  are  a  form  of  semiconductor  laser 
which  have  their  cavity  oriented  orthogonally  to  the  plane  of  the  wafer.  This  orientation 
necessitates  short  cavities,  requiring  highly  reflective  mirrors  and  a  relatively  high  gain/loss 
ratio  in  the  active  region.  The  resultant  superior  exit  beam  characteristics  (circular  beam 
output,  low  threshold),  the  capability  of  the  structure  to  be  arrayed,  and  the  tight  packing 
density  of  the  finished  lasers,  provide  strong  motivation  for  pursuing  the  growth  of  these 
structures.  This  thesis  details  the  design  of  an  optically  pumped  InGaAs  multiple  quantum 
well  periodic  gain  structure  VCSEL  with  a  950  nm  lasing  wavelength.  These  growths  are 
to  be  a  first  attempt  at  VCSEL  construction,  so  part  of  this  study  included  verification  of 
the  quality  of  the  back  mirror  growth  to  determine  if  it  met  the  design  specifications.  The 
research  included  construction  of  a  laboratory  configuration  to  measure  spectral  reflectivity 
and  optically  pump  the  finished  VCSELs.  This  provided  a  means  to  guide  the  growers  in 
the  fabrication  of  these  devices,  and  validate  la.sing  properties  in  the  finished  devices.  The 
VCSELs  were  characterized  by  spectral  reflectivity,  output  beam  polarization,  mode,  lasing 
wavelength,  and  optimal  pump  wavelength.  Analysis  of  the  characteristics  for  several 
VCSELs  obtained  from  the  University  of  Arizona,  and  the  back  mirror  grown  locally, 
illustrated  the  ability  to  use  measured  data  and  theoretical  spectral  reflectivity  calculations 
to  determine  the  quality  of  the  growths.  This  method  can  be  used  to  provide  a  guide 
for  subsequent  VCSEL  growths.  The  reflectivity  analysis  can  determine  if  the  growth 
parameters  have  been  roughly  achieved,  but  do  not  guarantee  the  quality  and  quantity  of 
lasers  on  the  wafer. 
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Design  and  Characterization  of  Optically  Pumped  Vertical  Cavity  Surface 

Emitting  Lasers 


I.  Introduction 

This  thesis  investigates  the  design  and  characterization  of  optically  pumped  vc-tical 
cavity  surface  emitting  lasers.  Research  into  the  design  of  Surface  Emitiing  Laser  (SEL) 
structures  has  increased  dramatically  over  the  past  decade.  Initially,  this  type  of  design 
was  explored  to  take  advantage  of  the  ability  to  array  the  SEL  in  two  dimensions  [29], 
allowing  tighter  packaging  of  individually-addressed  lasers.  Vertical  Cavity  SEL  (VCSEL) 
designs  have  been  aggressively  pursued  o^'er  the  past  few  years  because  VCSELs  have  ad¬ 
ditional  inherent  advantages.  The  VCSEL  design  exhibits  better  exit  beam  quality,  is  of 
smaller  size,  and  can  entail  less  complex  fabrication.  These  advantages,  coupled  with  the 
increasingly  widespread  access  to  the  necessary  fabrication  technology  to  make  VCSELs 
(e.g.  Molecular  Beam  Epitaxy  (MBE)  and  Metal-Organic  Chemical  Vapor  Deposition 
(MOCVD))  motivated  the  research  in  this  area  over  the  past  ten  years.  VCSELs  are  an 
emerging  technology,  and  may  provide  dramatic  improvements  in  surface  emitting  thresh¬ 
olds,  power  output,  and  optical  channel  capacity.  To  now,  commercial  applications  have 
driven  the  research.  Military  applications  in  higher  output  power,  beam  steering,  and 
heads-up  displays,  for  example,  may  proceed  more  briskly  with  the  development  of  an 
in-house  VCSEL  fabrication  and  test  capability.  This  research  provides  an  initial  design 
and  characterization  capability. 

1.1  Background 

The  typical  VCSEL  consists,  as  with  other  semiconductor  lasers,  of  a  semiconductor 
gain  medium  (III  -  V  compound)  and  a  set  of  mirrors.  The  design  of  both  the  gain  medium 
and  the  mirrors  is  critical  to  output  performance.  Although  it  is  an  integrated  structure, 
thp  VCSEL  can  be  thought  of  as  having  three  regions  (see  Figure  1.1).  The  top  and  bottom 
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regions  are  Distributed  Bragg  Reflector  (DBR)  mirrors,  and  the  middle  region  is  the  laser 
cavity. 
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Figure  1.1.  Vertical  Cavity  Surface  Emitting  Laser  (VCSEL  Geometry) 


1.1.1  DBR  Mirrors  The  advantage  of  MBE/MOCVD  lay  down  technology  is 
that  these  methods  each  can  deposit  the  mirrors  and  cavity  in  one  uninterrupted  process. 
This  is  done  by  forming  the  mirrors  of  multi-layer  films  (see  Figure  1.1  inset).  The  change 
of  refractive  index  between  layers  in  the  multi-layered  stack  results  in  a  reflected  and 
transmitted  component  of  the  incident  light  at  each  interface.  The  multiple  reflections  and 
transmissions  in  the  stack  can  be  engineered  to  selectively  transmit  or  reflect  at  desired 
frequency  bands  [26].  MacLeod  outlines  criteria  governing  the  reflectances  at  each  interface: 

•  The  intensity  reflectance  of  light  at  each  interface  is  given  by  [(1  -  p)  /  (1  -h  p)]^  where 
p  is  the  ratio  of  the  refractive  indices  of  the  media  on  each  side  of  the  interface. 

•  There  is  a  phase  shift  of  180®  when  the  reflectance  is  in  a  medium  with  a  lower 
refractive  index  than  the  adjoining  medium.  The  phase  shift  is  0®  for  higher-  to 
lower-index  transitions. 
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•  The  components  of  light  from  the  multiple  reflections  in  the  stack  constructively 
interfere  (add)  if  the  relative  phase  shift  is  a  multiple  of  360®  and  destructively 
interfere  if  the  relative  phase  shift  is  an  odd  multiple  of  180®. 

In  the  DBR,  high  reflectance  is  accomplished  by  making  the  layers  in  the  stack  pre¬ 
cisely  one-quarter  wavelength  (A/4)  thick  for  the  wavelength  to  be  reflected.  The  light 
incident  at  one  interface  will  traverse  this  thickness  and  partially  reflect  off  the  next  inter¬ 
face,  and  each  succeeding  interface.  The  layers  are  A/4  thick,  which  produces  a  tt  phcise 
shift  between  successive  reflections,  and  there  is  an  added  tt  phase  shift  at  the  lower- 
to  higher-index  interfaces;  as  these  components  return,  they  will  be  in-phase  and  con¬ 
structively  interfere,  producing  an  overall  high  reflectance  at  the  wavelength  the  stack  is 
designed  for.  Wavelengths  of  incident  light  which  are  significantly  longer  or  shorter  than 
this  wavelength  will  have  less  dramatic  reflectance  from  the  stack. 

The  DBR  mirror  design  produces  highly  reflective  mirrors.  We  must  consider  more 
than  just  the  reflectance,  however,  for  the  mirrors  to  form  a  good  semiconductor  laser 
mirror.  The  refractive  index,  absorption  coefficient,  and  lattice  matching  considerations 
all  impact  the  choice  of  mirror  materials.  These  factors  will  be  addressed  to  a  greater 
depth  in  later  sections. 

1.1.2  VCSEL  Cavity  The  cavity  for  the  VCSEL  contains  the  gain  medium,  where 
the  population  inversion  of  carriers  allows  stimulated  emission  to  produce  the  exiting  beam. 
This  cavity  can  be  completely  composed  of  active  material,  but  it  is  now  more  common 
to  see  a  cavity  which  uses  a  periodic  gain  structure.  That  is,  the  active  regions  of  the 
cavity  reside  only  at  the  places  in  the  cavity  where  the  electric  field  standing  wave  will 
form  a  maximum  (see  figure  1.2).  These  active  regions  can  be  formed  of  bulk  material  or  of 
multiple  quantum  wells,  with  barriers  interspersed,  to  take  advantage  of  quantum  effects. 
The  choice  of  materials  in  this  region  is  affected  by  the  desired  gain  for  the  region,  the 
absorption  of  the  spacer  and  the  active  regions,  the  lasing  wavelength,  and  lattice  matching 
considerations. 
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Figure  1.2.  Periodic  Gain  Structure  (PGS)  Cavity 


1.1.3  VCSEL  Characteristics 

Reflectivity  Together,  the  cavity  and  the  mirrors  form  an  integrated 
Fabry-Perot  etalon.  In  the  VCSEL,  the  Fabry-Perot  etalon  structure  provides  an  effective 
energy  storage  mechanism.  The  spectral  reflectivity  curve  is  characterized  by  regions  with 
high  transmittance  and  reflectance,  and  optically  pumped  VCSEL  structures  exploit  these 
areas  of  the  curve. 

A  typical  curve  for  a  bulk  GaAs  VCSEL  illustrates  the  Fabry-Perot  effect  in  Figure 
1.3.  The  resonance  at  880  nm  is  flanked  by  highly  reflective  regions  (~  99.9  %).  This  high 
reflectivity  is  achieved  with  less  than  25  pairs  of  AlAs/GaAs  A/4  layers  in  the  top  and 
bottom  mirror  stacks.  In  optical  pumping,  another  laser  is  used  to  generate  carriers  in 
the  VCSEL.  Thus,  the  mirrors  must  be  highly  transmissive  at  the  “pumping”  wavelength 
to  allow  the  laser  light  to  reach  the  gain  medium.  For  this  reason,  carriers  are  injected 
at  a  wavelength  (the  pump  wavelength)  near  the  edge  of  the  highly  reflective  stop  band, 
where  the  structure  is  highly  transmissive.  In  this  case,  a  likely  pump  wavelength  would 


be  near  820  nm.  The  curve  below  was  constructed  using  the  transmission  matrix  approach 
discussed  in  Chapter  II  and  the  FORTRAN  programs  included  in  Appendix  A. 
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Figure  1.3.  Fabry  Perot  Effect  in  Typical  VCSEL 

The  reflectivity  curve  illustrated  in  Figure  1.3  is  dependent  on  the  uniformity  of  the 
layers  in  the  mirrors,  the  cavity,  and  the  substrate,  and  the  index  of  refraction  of  the 
materials  in  these  regions.  The  ability  of  the  manufacturer  of  the  VCSEL  to  control  these 
parameters  will  determine  how  closely  the  flnished  product  meets  the  reflectivity  curve. 

Polarization  Strain  in  the  interfaces  between  layers  in  the  active  regions 
and  in  the  mirrors  can  affect  the  output  characteristics  of  the  VCSEL.  In  InGaAs  quantum 
well  active  regions,  strain  can  produce  polarized  output  [6].  Compressive  strain  tends  to 
produce  TE  polarization,  while  tensile  strain  favors  TM  polarization  [10].  For  active  regions 
constructed  of  AlxGai_,As  layers,  the  close  lattice  match  in  these  ternaries  over  the  range 
of  X  concentrations  provides  minimal  strain.  Polarization  measurements  for  bulk  GaAs 
active  region  VCSELs  indicate  no  favored  direction  in  the  linear  polarizations  [21]. 

Mode  Selection  Single  longitudinal  mode  operation  is  a  highly  valued 
inherent  feature  of  VCSELs,  even  though  typical  in-plane  semiconductor  lasers  can  have 
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natural  mode  spacings  on  the  order  of  hundreds  of  GHz.  The  small  cavity  size  (~  3/zm) 
makes  the  Fabry- Perot  structure  highly  mode  selective.  The  actual  mode  spacing  is  difficult 
to  calculate,  since  the  DBR  mirrors  have  penetration  depths  which  contribute  to  an  effective 
cavity  length  which  can  be  much  greater  than  the  designed  cavity  length  [20].  The  effective 
cavity  length,  however,  is  normally  small  enough  to  preclude  higher  order  longitudinal 
modes  in  VCSEL  structure. 

Higher  order  transverse  modes  have  been  noted  in  several  studies  [17].  Jewell  noted 
that  higher  pump  intensities  and  slightly  misaligned  input  beams  were  required  to  excite 
these  modes  [21].  The  transverse  mode  order  VCSELs  exhibit  depends  on  the  transverse  di¬ 
mensions  of  the  structure,  larger  transverse  dimensions  supporting  higher  order  transverse 
modes  [8]. 


Power  The  power  from  theses  devices  depends  primarily  on  the  cooling 
efficiency  and  the  power  efficiency  of  the  geometry.  Output  power  is  low  in  comparison  to 
other  types  of  semiconductor  lasers,  due  to  the  small  geometries  involved,  which  limit  the 
gain  volume  and  provide  heat  dissipation  problems  [19].  Pulsed  output  power  has  been 
reported  up  to  .5  W  and  continuous  wave  (CW)  power  above  9  mW  [9]. 

1.2  Problem  Statement 

This  thesis  examines  the  design  techniques  for  building  VCSELs,  and  the  character¬ 
ization  of  these  designs  via  surface  reflectivity,  output  polarization,  frequency,  and  beam 
quality. 

1.3  Scope 

The  scope  of  the  effort  was  to  design  optically  pumped  VCSELs,  and  develop  a 
laboratory  capability  to  guide  the  fabrication,  check  quality,  demonstrate  lasing  in  the 
final  product,  and  measure  completed  VCSEL  output  beam  characteristics  and  spectral 
reflectivity.  The  design  considerations  examined  for  this  thesis  can  be  applied  to  produce 
optically  or  electrically  pumped  surface  emitters.  This  thesis  purposely  focuses  on  the 
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optically  pumped  structures,  as  the  means  of  producing  VCSELs  for  this  effort  was  limited 
to  optically  pumped  devices. 

1.4  Approach/Methodology 

This  research  had  two  ends  to  achieve;  the  design  of  sample  vertical  cavity  surface 
emitters;  and  to  establish  a  means  to  characterize  them  via  direct  laboratory  measurement. 

The  design  part  involved  building  structures  to  be  optically  measured.  These  struc¬ 
tures  were,  basically,  prospective  parts  of  a  finished  VCSEL  design.  The  initial  effort 
was  segmented  to  provide  a  way  to  evaluate  the  growth  capabilities  of  the  available  MBE 
system  without  having  to  grow  entire  VCSELs. 

Most  current  designs  use  distributed  Bragg  reflector  mirrors  made  of  AUGai_*As, 
the  high  index  layers  using  a  low  x  concentration,  and  the  low  index  layers  using  a  high  x 
concentration.  AlAs/GaAs  is  a  conventional  choice  for  the  DBR  mirrors  for  InGaAs  lasers, 
and  was  chosen  because:  it  would  provide  a  predictable  means  of  determining  fabrication 
quality;  the  MBE  capability  available  to  us  was  ready  to  produce  this  combination  without 
modification;  and,  successful  VCSELs  have  been  grown  using  this  material  before  [17]. 

The  back  mirror  was  measured  to  see  if  it  exhibited  an  adequate  match  to  the  desired 
spectral  reflectivity  curve.  Once  the  back  mirror  design  was  complete,  the  full  VCSEL 
structure  was  submitted  for  fabrication.  The  gain  medium  used  a  periodic  gain  structure 
with  multiple  quantum  wells  at  the  peaks  of  the  standing  wave  electromagnetic  mode  in 
the  cavity. 

The  characterization  phase  involved  constructing  a  laboratory  bench  configuration 
to  provide  a  means  to  measure  reflectivity,  polarization,  lasing  wavelength,  and  beam  cross 
section.  The  laboratory  configuration  is  discussed  in  Chapter  IV. 

1.5  Sequence  of  Presentation 

Chapter  II  introduces  the  theory  used  to  design  the  structures  measured  for  this 
thesis.  The  VCSEL  and  back  mirror  designs  produced  by  this  effort  and  options  for 
electrical  stimulation  and  other  improvements  is  contained  in  Chapter  III.  Chapter  IV 
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provides  details  on  the  laboratory  configuration  used  to  measure  the  growth  structures. 
Chapter  V  details  the  results  of  the  testing,  and  Chapter  VI  summarizes  the  findings  and 
recommendations . 
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11.  Theory 


Building  a  vertical  cavity  surface  emitting  laser  follows  the  design  steps  for  larger 
conventional  lasers.  The  scale  and  orientation  of  the  VCSEL  cavity  produce  the  design 
challenge.  The  first  VCSELs  were  built  in  the  early  1970’s  [29].  The  initial  structures 
involved  three  distinct  regions,  the  top  and  bottom  mirrors  and  the  gain  region.  The 
application  of  MBE  and  MOCVD  to  the  process,  though,  has  removed  most  of  the  anxiety 
over  the  transition  between  the  cavity  and  the  mirrors.  The  result  is  a  one-step  growth 
which  produces  an  integrated  mirror/cavity  Fabry-Perot  etalon.  This  section  discusses  the 
initial  design  considerations  involved.  This  order  of  presentation  elaborates  on  the  design 
choices  used  in  building  the  InGaAs  quantum  well  VCSELs  for  this  effort,  but  is  readily 
adaptable  to  bulk  active  medium  designs. 

2.1  Efficient  Vertical  Cavity  Lasers 

Although  the  primary  concern  in  growing  the  samples  for  this  effort  was  not  efficiency, 
the  ability  to  incorporate  efficient  designs  in  early  VCSELs  will  benefit  later  designs.  Not 
all  applications  for  VCSELs  may  require  the  enhanced  efficiency  these  features  bring;  but 
a  review  of  them  provides  a  clearer  insight  to  the  design  options  available. 

The  goal  is  to  minimize  the  threshold  current  density,  ;  the  minimum  current 
density  required  for  lasing  to  occur.  J**  may  be  expressed  in  terms  of  the  nominal  current 
density,  the  active  region  thickness,  and  the  internal  quantum  efficiency.  The  nominal 
current  density,  inm ,  is  the  current  density  necessary  to  excite  a  1  /xm-thick  active  layer 
uniformly,  given  unity  quantum  efficiency  [30].  3,^  is,  then,  given  by  the  expression: 

(2.1) 

Vi 


where: 

Jth-  threshold  current  density  (per  unit  area) 
J„0m-  nominal  current  density  (per  unit  volume) 
d:  total  active  region  thickness 
Tii'.  internal  quantum  effidency 
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Applying  this  equation  to  the  VCSEL  geometry  will  optimize  threshold  current  den¬ 
sity,  and  provide  the  means  to  choose  the  correct  design  parameters  for  the  VCSEL. 

The  threshold  for  lasing  in  the  cavity  is  reached  when  the  light  traversing  the  cavity 
makes  a  round  trip  without  attenuation.  The  expression  for  this  condition  can  be  found 
in  most  introductory  semiconductor  laser  texts  [30:Chapter  7]: 


or,  taking  the  natural  logarithm  and  rearranging: 


(2.2) 


Tgth  =  a  -I-  —  In 


RkR 


(2.3) 


where: 

reflectivity  of  the  back  mirror 
Rj’.  reflectivity  of  the  front  mirror 
F:  confinement  factor 
gth’’  threshold  gain  of  the  region 
a:  loss  per  unit  length  of  the  region 
1:  length  of  the  cavity 

Equation  2.3  is  the  focus  of  this  chapter.  Dramatic  improvements  in  VCSEL  thresh¬ 
old  current  density  stem  from  relatively  straightforward  improvements  in  cavity  design, 
active  region  design,  and  mirror  selection.  All  of  these  improvements  reduce  to  an  action 
on  one  of  the  parameters  in  equation  2.3.  The  following  sections  discuss  the  underlying 
theory  behind  these  improvements  and  how  the  designs  for  the  structures  grown  for  this 
thesis  effort  benefit  from  them.  Recent  improvements  in  cavity  design  include  the  Periodic 
Gain  Structure  (PGS),  discussed  briefly  in  chapter  1.  Active  region  improvements  include 
Quantum  Well  (QW)  active  regions,  and  strained-layer  effects.  In  mirror  design,  the  use 
of  distributed  Bragg  reflector  mirrors  was  an  early  improvement  to  the  first  [29]  VCSEL 
designs,  but  a  description  of  the  theory  is  included  here  to  complete  the  treatment  of  the 
entire  VCSEL  structure. 
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2.2  Cavity  Design 

Ornate  cavity  designs  for  VCSELs  are  not  essential  to  lasing.  A  simple  block  of  bulk 
active  material  (e.g.  GaAs)  has  been  successfully  used  in  VCSEL  design  from  the  VCSEL’s 
inception.  This  type  of  cavity  design  is  now  referred  to  as  the  Uniform  Gain  Structure 
(UGS)  [11].  Periodic  Gain  Structures  (PGS)  have  become  more  popular  recently,  as  they 
can  theoretically  provide  up  to  a  factor  of  two  improvement  in  the  threshold  current  density 
without  changing  the  component  materials  of  the  VCSEL.  The  PGS  takes  advantage  of 
the  physical  dimensions  of  the  cavity  to  concentrate  the  active  region  at  the  peaks  of  the 
cavity  electromagnetic  mode  (see  Figure  1.2).  To  see  how,  we  need  to  look  more  closely  at 
the  factors  involved  in  equation  2.3. 

2.2.1  P:  Confinement  Factor  The  confinement  factor,  P,  determines  how  much 
of  the  electromagnetic  field  in  the  cavity  will  coincide  with  the  active  region.  The  confine¬ 
ment  factor  may  be  defined  in  three  dimensions  for  the  vertical  cavity  structure;  we  are 
interested  in  the  P  for  the  longitudinal  direction  of  the  cavity.  The  longitudinal  confine¬ 
ment  factor,  Pp,  is  the  ratio  of  the  magnitude  of  the  electric  field  in  the  active  region  to 
that  in  the  overall  cavity: 


r  _  Lcuve  E\z)dz 
UE^{z)dz 


(2.4) 


To  determine  Pp,  we  use  an  approximation  of  the  standing  wave  pattern  of  the  electric 
field  in  the  cavity  [11].  For  high  reflectance  mirrors  (which  is  usually  the  case  in  VCSELs): 


E{z)  =  Eocos{kz)  (2.5) 

For  the  PGS,  assuming  the  active  regions  are  placed  at  each  of  the  peaks  of  the 
electromagnetic  mode  in  the  cavity,  the  spacing  of  the  active  regions  is  A/2  (see  Figure 
2.1).  Defining  the  thickness  of  each  active  region  in  the  PGS  as  f,  and  the  totzd  active 
region  thickness  as  d  (which  is  the  sum  of  the  active  regions  t  in  the  cavity),  the  ratio  of 
d/t  is  the  number  of  periods  in  the  standing  wave  pattern.  This  ratio  is  the  same  as  the 
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ratio  of  the  cavity  length  to  one-half  the  lasing  wavelength  (j^).  Inserting  equation  2.5 
into  equation  2.4  and  integrating  over  the  thickness  for  one  period  [11]: 


rp 

rp 


ff,Eicos\kz)dz 

A72/A/2-^oCOs2(l:2)dz 


(2.6) 


Top 


Mirror 


Figure  2.1.  Periodic  Gain  Structure  (PGS)  VCSEL  Cavity  Geometry 


Corzine  defines  Pr  in  equation  2.6  as  the  relative  confinement  factor  [11].  For  the 
situation  described  above,  where  the  active  regions  lie  at  every  peak  of  the  standing  wave 
E-field,  we  can  equate  the  ratios  d/i  and 


d  I  d  t 

-  =  -7-  =>  j  =  -rjo  active  region  at  every  peak 
t  njzi  i  A/2 


This  provides  a  simple  equation  for  the  relative  confinement  factor  by  substituting  into 
equation  2.6. 
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(2.7) 


_  ,  ,  sin(Tf) 
(Trf) 


relative  confinement  factor 
total  active  region  thickness 
cavity  length 


The  factor  d/l  is  very  useful  in  the  design  of  the  cavity,  d/l  is  called  the  fill  factor  for 
the  cavity  and  represents  the  ratio  of  the  active  region  thickness  to  total  cavity  thickness. 
Note  that  from  equation  2.7  the  relative  confinement  factor  is  maximized  for  small  values 
of  d/l.  For  small  argument  form,  the  term  sinx/i  ~  1,  and  the  maximum  value  for  Fr  then, 
is  2  (we’ll  recall  this  later).  Now,  substituting  Frd//  from  equation  2.6  for  the  longitudinal 
confinement  factor  (F  =  Fp)  in  equation  2.3: 

^  d  1.  ,1, 

Tri/.hy  =  a  +  yln[-] 

^V,g,Hd  =  a/  +  ln[i]  (2.8) 

Where  R  is  the  geometric  mean  reflectivity  of  the  front  and  back  mirrors 
This  equation  provides  the  relationship  between  the  gain  of  the  PGS  active  medium  and 
the  losses  in  the  cavity  and  mirrors  for  threshold  gain. 

2.2.2  Cavity  losses  The  losses  can  be  further  subdivided  into  the  specific  loss 
factors  for  each  region  of  the  cavity  the  light  passes  through  as  it  makes  the  round  trip 
between  reflections.  We  can  subdivide  the  loss  term,  al  into  the  loss  in  the  cavity  itself 
and  the  loss  for  the  light  penetrating  the  mirror  layers  adjacent  to  the  cavity  (see  Figure 
2.1)  [9]. 


^r9thd  —  Ojl  +  Omlm  "b 


2-5 


where: 


d:  total  active  region  thickness 
lm‘  penetration  depth  of  mirrors 
TrSud:  active  region  gain 
ail:  loss  in  the  cavity 
amlm-  mirror  penetration  loss 

The  break-out  of  the  a, 7  loss  term  in  equation  2.9  includes  the  loss  in  the  active  and 
passive  regions  of  the  cavity,  and  loss  from  scattering  at  the  sidewalls  of  the  VCSEL  [9]: 


Og  —  "b  (l  Ep)Orgp  -j-  Ofgc 


(2.10) 


with: 

Ogg,:  the  loss  per  unit  length  in  the  active  region 
aipi  the  loss  per  unit  length  in  the  passive  region 
a,c:  scattering  loss  for  the  cavity 

The  free  carrier  loss  for  GaAs  in  the  active  elements  of  the  structure,  is  given  by 
Casey  [7:pl75]  as: 


Oia  «  3  X  10-‘®n-i-7  X  10-^®p  (2.11) 

For  n  «  p,  as  in  un-doped  GaAs  this  reduces  to  Oia/n  =  10“^^. 

Loss  mechanisms  are  less  well-documented  for  InGaAs.  Coldren  suggests  loss  factors 
can  be  as  high  as  100  cm“*  for  large  carrier  densities.  However,  he  later  cites  a  typical 
active  region  loss  for  InGaAs  as  30  cm“^  [9:S112). 

Passive  region  loss  and  scattering  loss  are  device  dependent,  for  Al*Gai_jAs  passive 
regions,  various  authors  have  used  from  1  to  10  cm“‘  [28]  [9].  We  will  use  10  cm~^  when 
we  require  a  loss  mechanism  for  later  calculations.  Scattering  losses  are  also  sparsely 
documented.  Corzine  uses  an  a,e  =  20cm“*  as  typical  [11]. 

2.2.S  Ja  Design  Equations  From  equation  2.1  and  equation  2.9  Corzine  develops 
a  useful  design  expression  for  the  threshold  current  density.  To  do  so  we  solve  equation 
2.9  for  d  and  substitute  this  expression  into  equation  2.1. 
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_  *^nom^  _  *^nom  “I"  Otfn^m  “f"  1 

m  Tji  TrQth 

_  ^nom  (^m^m  “H  ^^('fl])  ^m^m  'r 

Vi  ^rQth  ^^[^]) 


(2.12) 


By  rearranging  and  using  the  expression  for  Fp  (equation  2.3)  the  resultant  threshold 
current  design  equation  falls  out.  From  equations  2.3  and  2.9: 


Thus; 


^p9th^  —  “I"  “h 


~h  Qm^m  d*  l^[fi])  _  Fp^tfc/  _  ^pdth 
(^mffn  "h  ^^[^])  ^p9th^  f'p^tA 


Substituting  this  into  equation  2.12: 


_  Jnom  1  Ml/  )  Tpgth 

5**  aJl+L(i/R)  ~ 

■/nom  ln(l/iZ)  1  L  Oj  1-1 

9th  VmVi  Tr  ^p9th' 


(2.13) 


Equation  2.13  separates  the  components  of  the  threshold  current  density.  It  is  useful 
to  examine  the  components  of  equation  2.13  in  more  detail. 
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^  normalization  constant  Jo(R-)  and  represents  the  current 

density  which  will  normalize  the  threshold  current  density  Jth  for  no 
internal  losses  (a,  =  0)  and  for  a  relative  confinement  factor  F,  =  1 
(Note:  this  infers  the  fill  factor  d/l  =  1,  corresponding  to  Uniform 
Gain  Structure  (UGS)  active  medium) 

Fr  The  relative  confinement  factor,  as  noted  before,  has  a  maximum 
value  of  2,  indicating  a  potential  for  F,  to  halve  the  threshold  current 
density. 

ratio  of  the  nominal  current  density  to  the  threshold  gain.  By  mini¬ 
mizing  this  factor,  we  may  minimize 

Efficiencies:  Although  it  is  unnecessary  to  separate  the  efficiencies  of  the  cav¬ 
ity,  mirror,  and  active  regions,  it  can  be  helpful  in  relating  design 
efforts,  to  compare  mirror  efficiencies  for  a  design  using  varying  con¬ 
centrations  of  Aluminum,  for  example.  The  various  efficiencies  for 
the  structure  are  explained  below 

rim-  the  conversion  efficiency  of  the  mirrors  is  given  by  the  equation: 

ln(l/i?) 

a,n/„-bln(l/il) 

and  is  1  for  no  mirror  penetration  loss  (om^m  =  0)- 

Tii'.  internal  quantum  efficiency  (for  the  active  region  only) 

Tic’,  the  cavity  efficiency,  together  with  the  mirror  efficiency  T)m  pro¬ 
vide  the  external  efficiency  of  the  structure  (ratio  of  output  loss 
though  mirrors  to  total  loss  (cavity  and  mirrors)): 


Vtx  = 


“b  ^rrJ-m  "b 


"b 


(2.14) 
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2.3  Active  Region  Design 

VCSEL  active  region  design  improvements  have  fed  off  improvements  in  active  regions 
for  in-plane  semiconductor  lasers.  The  addition  of  multiple  quantum  wells  and  strained 
layers  led  to  improved  threshold  current  densities  in  addition  to  those  attainable  by  the 
PGS  structure  presented  in  the  previous  section.  It  is  instructive  to  first  examine  the 
contribution  of  plain  bulk  active  regions  to  the  PGS  structure,  then  note  how  QWs  can 
improve  performance: 

2.3.1  Bulk  Active  Regions  For  bulk  active  regions,  to  obtain  the  minimum  J,/, 
we  need  to  determine  the  minimum  normalization  constant  Jo(R).  To  this  end,  we  need 
to  maximize  ■qmi  »?•»  and  R,  and  minimize  J„om/g«/>  in  equation  2.13.  Of  these  parameters, 
active  region  design  influences  i„omlgth  most. 

To  illustrate,  the  process  for  a  GaAs  active  region  follows.  A  derivation  for  the  linear 
relationship  between  gain  and  carrier  density  for  GaAs  can  be  found  in  Yariv  [33:p564]. 

g  =  a{n-ntr)  (2.15) 


with: 

n  carrier  density  in  the  material 
Utr  carrier  density  necessary  for  transparency 
a  curve  fitting  constant 


From  Mcllroy  [27]  for  small  dopant  levels  (n  w  p)  the  nominal  current  density  can 


be  given  by  the  equation: 


(2.16) 


with: 

e  electron  charge 
Bej j  recombination  rate 

By  solving  equations  2.16  and  2.15  for  n,  and  equating  the  resultant  expressions,  we 
can  obtain  an  expression  for  g  in  terms  of  J„om : 
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•^nom  —  ^  — 

g  =  a(n  —  ntr)  =>•  n  =  -  +  n^r 

a 

=•  S  =  “1 -  "-1  (21’') 

Corzine  assigns  the  following  values  to  the  constants  in  this  equation  [11]: 

a:  2.5  x  10“^®CTn^ 

Be}/'-  lQ~^°cm^/s 
ritr'-  1.5  X  10'*cm“^ 

With  these  assignments,  a  graph  of  the  gain  vs  values  of  J„^  illustrates  their  inter¬ 
dependence  (see  Figure  2.2) 


Figure  2.2.  Gain  vs  Nominal  Current  Density  Curve  for  GaAs 


By  dividing  both  sides  of  equation  2.17  by  J„(,m  we  can  also  plot  g/J„om  vs  Jnom  (see 
Figure  2.3) 

If  we  maximize  g/Jnom  from  figure  2.3,  this  fixes  Jnom')  we  can  then  extract  the 
corresponding  gain  from  figure  2.2.  As  a  first  estimate,  Corzine  suggests  the  value  for  the 
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Jfiom(10Ag  A/cin3) 

Figure  2.3.  g/Jnom  vs  J„om  Curve  for  GaAs 


gain  should  be  chosen  on  the  high  side  of  the  peak  in  figure  2.3  [11).  By  doing  so  the  gain 
will  help  to  minimize  the  effects  of  internal  loss  in  the  cavity  efficiency  (see  equation  2.14). 

d/1;  Fill  Factor  Choice  The  fill  factor  choice  is  the  first  point  in  the 
analysis  where  the  improvement  in  efficiency  is  dramatically  evident.  Using  equation  2.13, 
we  can  write  the  ratio  of  the  threshold  current  density  to  the  normalization  constant  as: 


Jth  _ 

MR)  ~  Tr 


(2.18) 


For  example,  in  GaAs,  we  can  determine  an  optimal  fill  factor.  Assume  reasonable 
internal  losses  (  Q,p  =  lOcm”*,  a,c  =  20cm~*:  see  Cavity  Losses  above).  Determine 
the  optimal  gain  by  choosing  a  Jnom  in  figure  2.3  which  produces  a  maximum  in  g/J„om- 
Now  use  this  J„ot„  to  determine  the  gain  from  figure  2.2.  The  result  is  a  J„om  between 
1.2  X  lO^A/cm^  and  1.8  x  UfAjcm^,  and  corresponds  to  a  threshold  gain  in  the  vicinity  of 
350  to  500  cm“*.  As  noted  before,  Corzine  [11]  advises  the  gain  should  be  chosen  near  the 
high  end  of  this  range,  to  minimize  losses.  Assigning  a  gain  of  500  cm~',  we  can  determine 
the  carrier  density  n  from  equation  2.15.  With  the  constants  a  and  n,r  as  defined  before. 
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the  gain  of  500  cm“^  implies  n  =  3.5  x  10‘®.  Now  insert  this  value  (assume  n  «  p  true  for 
un-doped  material)  into  the  expression  for  the  free  ci  jrier  absorption  in  equation  2.11: 

Oia  =  n  X  10”^^ 

=  35cm~^ 


Returning  to  equation  2.18,  if  we  plug  in  the  expressions  for  Oj  from  equation  2.10, 
and  Fp  from  equation  2.6: 


JtH 

MR) 


rpQia  -b  (1  -  rp)a,p  +  a, 

r,(d//)Q,a  +  (i-r,(d//))a.p  +  Q. 


Tr(d/l)ga 


Finally,  using  F,  from  equation  2.7: 


(2.19) 


Tr 


=  1  + 


sin(irj) 


We  can  plot  equation  2.19  vs  fill  factor  d/l. 

The  plot  (see  figuie  2.4)  illustrates  how  the  correct  choice  of  the  fill  factor  can  improve 
the  threshold  current  density.  Unit  gain  in  the  graph  represents  the  gain  for  a  UGS  GaAs 
active  region  with  no  loss.  The  dip  in  the  curve  represents  a  40  %  reduction  in  the  threshold 
current  density  for  a  periodic  Gain  Structure  (PGS)  including  loss  over  a  Uniform  Gain 
Structure  (UGS)  with  no  loss. 


2.3.2  Quantum  Well  Active  Regions  Quantum  well  active  regions  can  provide 
additional  reductions  in  threshold  current  density.  For  both  InxGai_xAs  and  GaAs,  the 
optimal  gain  of  the  QW  is  dependent  on  the  well  width.  A  brief  background  in  quantum 
well  theory  is  in  order  to  determine  how  to  make  informed  choices  concerning  quantum 
well  parameters. 
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Figure  2.4.  JthlJo(R)  vs  Fill  Factor(d//)  Curve  for  GaAs 

Transition  Energy  The  transition  energy  in  quantum  wells  is  comprised 
of  four  primary  contributors  [32]  (see  figure  2.5). 


E  =:  hu  —  Ea  +  -C'on/  +  RconJ  ~  Rtxc 


(2.20) 


where: 

h:  Boltzman’s  constant 
v:  frequency  of  the  transition 
Ea-  energy  gap  of  the  bulk  semiconductor 
Reonf-  confinement  energy  of  the  conduction  band 
^eonf-  confinement  energy  of  the  valence  band 
Etce-  exciton  absorption  energy 

Gap  Energy  The  energy  gap  for  In*Gai_jpAs  and  Al*Gai_,As  bulk 
active  regions  depend  on  the  concentrations  of  the  Indium  (or  Aluminum)  in  the  ternaries. 
The  In,Gai_,As  parameters  for  this  effort  were  obtained  from  Kolbas  [22],  whiie  the 
Al,Gai_eAs  parameters  were  gotten  from  Adachi  [2] 
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Figure  2.5.  Quantum  Well  Energy  Diagram 


Electron  Confinement  Energy  The  electron  confinement  energy  can 
be  calculated  using  the  infinite  well  approximation.  For  initial  calculations  we  take  the 
energy  difference  AiJe  between  the  conduction  bands  of  the  well  and  the  barriers  to  be 
infinite.  Modelling  the  energy  in  the  system  with  the  Schrodinger  equation;  we  have,  for 
the  conduction  band: 


V{z)u{r)- 


-^V=o(r)  =  Eu(r) 
2mJ  '  ^ 


(2.21) 


where: 

V(z):  potential  energy  function  confining  the  electrons  in  the  z  direction 
u(r):  wavefunction  for  the  electrons  in  the  well 
E:  energy  of  the  electron 
m* :  effective  mass  of  the  electron 

Applying  separation  of  variables  to  solve  the  equation,  the  eigenfunction  u(r)  can  be 
written  as  the  product: 


«(r)  =  V’t(ri)«(z)  (2.22) 

. .  .with  fj.  the  component  of  r  perpendicular  to  the  z  direction  (in  the  plane  of  the 
well).  By  inserting  this  relation  into  equation  2.21  and  extracting  the  z-dependent  portion 
we  obtain  a  1-D  wave  equation: 
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=»  V{z)  - 


2m! 


V{z)MrL)uiz)  -  — VVi(rx)u(0) 

=>  ^^(2)- ^V2]^t(rxM2) 

^^(^*(rx)) 


r  1  S^ujz))  1 
u{z)  6z^  V’t(rx) 


6x^ 

1  S^iiptirj.)) 


+ 


^t(rx) 

Now,  extracting  the  z-dependent  component: 


6yi 


Eipk{rA.)u{z) 

Etl}t{Tx)u{z) 


E 


(2.23) 

For  well  width  L,,  u{z)  must  vanish  at  the  walls  {z  =  ±1^2)  for  the  equation  to 
remain  bounded,  since 

lim  V{z)  —*  00 

with  this  boundary  condition: 


fcoslj^z  /=  1,3,5,... 
\sin/f;2  /=  2,4,6,... 


Solving  for  the  constant  E,  when  y(2)  =  0  (in  the  well) 


(2.24) 


E.  =  l^ 


2m:Ll 


I  =  1,2,3. 


(2.25) 


The  energy  for  the  x  and  y  components  in  equation  2.21  are  determined  by  solving 
for  V’(*’x)  ^  ^  two  dimensional  Bloch  wavefunction  [33:Chapter  15].  The  wavefunction  of 
the  electron  in  the  conduction  band  is  characterized  by  the  wave  vector  kx  and  a  Bloch 
wavefunction: 

V’(rx)  =  Ukx(rx)e^’'‘-*-  (2.26) 
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The  expression  of  the  overall  energy  of  the  electron  characterized  by  the  vector  kj.  as 
measured  from  the  bottom  of  the  band  is  [33]; 

£(k.)=0  (2.27) 

By  adding  the  energies  in  equation  2.25  and  equation  2.27,  the  overall  energy  above 
the  conduction  band  becomes: 


E:^j(kj^J)  =  E,  +  Eikj.)  = 


2m-  2m:Z2 


(2.28) 


Hole  Confinement  Energy  The  hole  quantization  problem  is  more 
difficult,  a  result  of  the  degeneracy  of  the  valence  subbands.  The  more  rigorous  methods 
used  to  characterize  the  energies  in  the  valence  band  for  light  and  heavy  hole  masses 
use  the  Kane  model  [32].  This  model  describes  the  valence  bands  using  basis  functions 
with  angular  momentum  symmetry  (J  =  ±§,±5).  This  equates  to  4-fold  degeneracy  at 
it  =  0.  The  basis  functions  are  described  by  the  Luttinger-Hamiltonian,  and  an  explicit 
development  is  contained  in  Weisbuch.  The  resultant  energies  for  the  light  and  heavy  hole 
confinement  energies  are  [32]: 


(2.29) 

(2.30) 

where: 

E^Vonj-  confinement  energy 

^cmj-  heavy  hole  confinement  energy 
7i,2:  Luttinger  parameters 
mo:  electron  rest  mass 


t2  2 

Ec^fik^.,1)  =  +  72)  +  t^(7i  +  72) 


0(kx,/)  = 


2moL^ 
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The  Luttinger  parameters  can  be  used  to  obtain  the  effective  mass  for  the  valence 
band  light  (m/^)  and  heavy  {nihh)  holes  in  the  strained  quantum  well  [12]: 


1 


m^A 

=  mo - 

'll  -72 

1 

m/A 

=  mo - 

7i  +72 

These  forms  are  convenient  if  the  Luttinger  parameters  can  be  obtained  for  the 
active  region  materials  involved,  one  source  for  Luttinger  parameters  in  III-V  compounds 
is  Landolt-Bornstein  [24]. 

Exciton  Energy  The  electron  and  hole  pairs  in  the  conduction  and 
valence  bands  are  correlated  in  their  motion  [32].  This  correlation  can  be  modeled  as  a 
simple  Coulomb  attraction  of  the  electron-hole  pair  (also  referred  to  as  an  exciton).  The 
effect  of  this  interaction  produces  the  exciton  binding  energy  E^xe  equation  2.20.  The 
exciton  effect  is  not  significant  in  lasers,  where  high  injection  of  charged  carriers  screens 
the  Coulomb  interaction  [5].  For  this  reason  we  will  assume  Etxc  ~  0. 

The  total  transition  energy  is  the  sum  of  the  energies  in  equation  2.20.  With  E^xc  ~  0: 


hu 


Eg  -f- 


2m,  V 


LV 


(2.31) 


where: 

Eg-  energy  gap  of  the  bulk  well  material 
L,:  well  width 

m,:  reduced  mass  of  the  electron  -  hole  pair 

1:  1,2,... 


Quantum  Well  Barriers  The  barriers  must  have  a  sufficiently  large  en¬ 
ergy  gap  above  that  of  the  wells  to  contain  the  carriers.  In  addition  the  barriers  must 
separate  the  wells  to  prevent  inter- well  coupling,  which  can  negate  quantum  effects  [32]. 
Jewell  [19]  suggests  barrier  widths  of  100  A,  and  Corzine  [10]  recommends  the  conduction 
band  barrier  height  be  kept  100  meV  above  the  well. 

Strain  Effects  Lattice  strain  in  the  active  layers  can  alter  the  dispersion 
of  the  valence  bsind  energies.  For  GaAs  and  In^Gai.^As  active  regions  on  GaAs  substrates, 
compressive  strain  (see  figure  2.6  [10])  tends  to  elevate  the  heavy-hole  valence  band  above 
the  light-hole  valence  band.  The  result  is  less  band  mixing  in  the  valence  band.  This 
makes  the  bands  more  parabolic  and  decreases  the  density  of  states  in  the  upper  valence 
bands.  The  lightened  upper  valence  band  density  of  states  reduces  the  current  density  for 
transparency  (Jtr),  lowering  the  threshold  current  density. 


WUMDliMd 


DOS:  Dmiiiy  of  SlaiM 
Hx:  Ibny  Hde  Sob-baMa 
LKUgteHoloSob-bttds 


Figure  2.6.  Compressive  Strain  Effects  on  Valence  Sub-Bands  (GaAs)  [After  Corzine[10]] 


Strain  for  InGaAs  weUs  on  GaAs  substrate  is  compressive  for  any  Indium  z  con¬ 
centration.  In.Gai.^As  wells  on  InP  substrate  produce  either  tensile  (for  x  <  .532)  or 
compressive  (z  >  .5321)  strain. 
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Critical  Thickness  It  is  important  in  strained  layer  design  to  be 
certain  the  layer  thickness  does  not  exceed  the  critical  layer  thickness.  The  critical  layer 
thickness  for  the  lattice  is  a  function  of  the  strain  and  the  layer  thickness.  Exceeding  this 
thickness  forms  dislocations  in  the  growth  surface  which  propagate  to  the  interface,  forming 
uneven  layer  growth  [22].  If  using  In*Gai_j;As,  to  achieve  the  desired  transition  energy,  the 
designer  must  pick  a  well  width  L,  and  an  Indium  concentration  x.  The  x  concentration 
affects  Eg,  and  m*  in  equation  2.31.  It  also  affects  the  strain  between  layers,  so  given 
a  well  width,  there  will  be  a  limit  on  x  concentrations  which  will  produce  layers  below  the 
critical  layer  thickness.  Fritz  characterizes  In^Gai_xAs  for  critical  thickness  vs  strain  for 
strain  values  of  0.5  %  to  3  %  and  for  layer  thicknesses  from  50  to  500  A[15]. 

Polarization  Polarization  of  the  output  beam  is  usually  attributed 
to  strain  in  the  active  regions  [14].  The  strain  in  InGaAs/GaAs  and  InGaAs/InGaAsP 
quantum  wells  has  been  shown  to  produce  significant  polarization  effects.  Transverse 
electric  (TE)  peak  gain  is  usually  defined  as  corresponding  to  the  gain  seen  by  an  optical 
field  with  its  electric  field  polarized  in  the  plane  of  the  well.  Transverse  magnetic  (TM) 
peak  gain  corresponds  to  the  electric  field  perpendicular  to  the  plane  of  the  well.  Ahn 
and  Chuang  [3]  link  TM  gain  to  the  relative  dominance  of  light  hole  sub-bands  in  the 
quantum  well;  while  they  attribute  TE  gain  to  heavy  hole  sub-band  dominance.  For  this 
reason,  compressive  strain  as  shown  in  figure  2.6  tends  to  favor  TE  gain.  Tensile  strain 
has  been  shown  to  raise  the  light  hole  sub-bands  in  InGaAs,  and  produced  predominantly 
TM  polarized  light  [10]. 

2.4  Mirror  Design 

The  distributed  Bragg  reflector  mirror  design  has  been  prevalent  in  VCSEL  struc¬ 
tures  since  the  late  70’s.  This  section  describes  the  transmission  matrix  approach  to 
calculating  mirror  reflectivities,  and  discusses  the  effects  of  the  penetration  depth  of  the 
electromagnetic  mode  into  the  cavity. 

2.4- 1  Transmission  Matrix  The  transmission  matrix  approach  to  reflectivity  cal¬ 
culations  models  the  multiple  reflections  in  the  multilayer  medium  as  two  lumped  optical 
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fields:  one  traveling  in  the  positive  direction,  and  one  in  the  negative,  or  reflected,  direc¬ 
tion.  The  positive  going  electric  and  magnetic  field  vectors  in  the  m"*  layer  are  denoted 
by  E+  and  H+.  The  reflected  fields  are  denoted  by  E“  and  H~ .  The  total  field  in  the 
m‘*  layer  can  then  be  written  [25]: 

Em  =  E;],  -h  E~ 

H„  =  /Xm(E+-E-)  (2.32) 

Where  the  final  equality  uses  the  (possibly  complex)  effective  refractive  index  {^lm) 
of  the  m**  layer  to  relate  the  E  and  H  fields.  The  effective  indices  in  the  mirror  layers  is 
derived  from  the  complex  refractive  index  (n^)  and  the  angle  of  incidence  (^m): 


=  ^m/cos^n,,  for  the  p-component  of  polarization 
=  Tim  cos  for  the  s-component  of  polarization  (2.33) 

6m  can  of  course  be  related  to  the  external  angle  of  incidence  {6o)  (see  figure  2.7) 
and  the  angle  of  the  beam  in  the  substrate  (^,„4)  by  Snell’s  law.  The  phase  of  the  light  as 
it  traverses  the  m*'*  layer  is  determined  by  the  equation: 

C  _  2  TT  Tim  dm  cos  ^m  fn 

dm  —  ^  (2.34j 

where: 

6m'  phase  of  the  field 
dm:  physical  thickness  of  the  layer 
A:  free  space  wavelength  of  the  light 

Following  the  matrix  method  developed  by  Heavens,  [18],  we  apply  boundary  con¬ 
ditions  to  the  E  and  H  fields  incident  at  the  interfaces.  Namely,  these  fields  must  be 
continuous  across  the  boundaries.  The  matrix  solution  to  this  boundary  value  problem  is 
of  the  form: 
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incktont  bMm  MtgctM  Compoo«nt 


Figure  2.7.  Geometry  for  Reflected  Fields  in  Multilayer  Stack 


cos  6m  sin  6m 
iHm  sin  6m  cos6m 


In  this  form,  we  can  easily  extend  the  application  to  an  N  layer  stack; 


(2.35) 


(2.36) 


with  Eo  and  Ho  the  incident  electric  and  magnetic  fields  and  the  characteristic 
matrix: 


cos  6m  ^  sin  6m 
tfim  sin  6m  COs6m 

We  are  interested,  for  the  VCSEL  mirrors,  in  the  reflectivity  in  a  stack  composed 
of  quarter- wave  thicknesses  (each  layer  has  dm  =  ^/^)-  For  this  thickness  the  phase  6m 
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becomes: 


(2.37) 


27rn„d„,cos0„  ircos^^ 
j - -  —5— 

For  normal  incidence,  {6^  =  0))  this  results  in  a  value  of  =  f  •  So  the  characteristic 
matrix  Mm  reduces  to: 


Mm  = 


0  -i- 


normal  incidence,  dm  =  A/4 


From  this  matrix  expression,  we  can  determine  the  exiting  fields  from  the  quarter- 
wave  stack;  for  N  layers  and  the  substrate: 


/ 

\ 


Eo 

Ho 


N+Mub 

n 


m=l 


E,ub 

H$ub 


N+sub 

=  n 


m=l 


0 


l>m  0 


1 

M»ui 


E. 


sub 


(2.38) 


The  reflectivity  is  then  given  by  [25]: 


fioEo  —  Hq 

2 

f^o  ~  Hof  Eo 

fioEo  +  Hq 

/Xo  +  Ho!  Eo 

(2.39) 


Once  the  product  of  the  matrices  in  equation  2.38  is  determined,  the  ratio  HqIEq  = 
y,  where  Y  is  the  optical  admittance  follows.  This  results  in  two  ciises  (for  N  odd  or 
even): 


y  = 


1  Q 

N  odd 

- — 173 - 

N  even 

■  ft 


And  the  reflectivity  is  a  straightforward  calculation  [25]: 


/XQ-y  ^ 

/xo  +  y 


(2.40) 


(2.41) 
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The  choice  of  materials  in  the  layers  dictates  the  values  of  Y.  The  object  is  to  have 
Y  »  (iQ  OT  Y  «  hq  to  make  J?  cs  1;  to  do  so  one  must  choose  the  layers  adjacent  to 
the  substrate,  and  incident  medium,  so  that  the  change  in  refractive  index  is  majdmized 
at  these  interfaces.  This  will  determine  whether  N  will  be  even  or  odd  in  equation  2.40. 
The  number  of  even  or  odd  layers  will  then  determine  the  admittance  and  reflectivity  to 
be  attained. 

2.4-2  Refractive  Indices  The  VCSEL  mirrors  in  this  thesis  effort  used  varying 
concentrations  of  Alj,Gai_*As.  The  material  parameters  based  on  x  concentration  for 
Al*Gai_jpAs  were  obtained  from  Adachi  [2]. 

The  refractive  index  is  taken  as  the  square  root  of  the  real  part  of  the  dielectric 
constant: 


n(w)  ~ 


with: 

n(u)):  frequency  dependent  index  of  refraction 

Ci(a;):  frequency  dependent  real  part  of  dielectric  constant 

And  ei(w)  for  AlGaAs  can  be  expressed  as  [2] 


ex{u)  =  Ao(a:)  |/(x)  +  ^  [Eo/iEo  +  Ao)f^  /(x.o)}  +  ^o(x)  (2.42) 

In  equation  2.42  the  components  are  defined: 


fix)  =  X  ^  [2-(l  +  x)*^^-(l-x)^^^] 

X  =  hu/ Eq 

X,o  =  ftu;/(Eo  +  ^o) 

where: 

Eq:  direct  band  gap  energy 
Aq:  spin  orbit  split-off  energy 
Ao,Bo:  x-dependent  constants 
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These  parameters  may  be  expressed  in  terms  of  the  x  concentration  of  the  material 


as  [2]: 


Eo  =  1.425+  1.155a: +  0.371^  eV 
fJo  +  Ao  =  1.765+  1.115X  + 0.37x2  eV 
Aq  =  6.30  +  19.0x 


Bo  =  9.4  -  10.2x 


2.4-3  Reflectivity  Calculations  The  above  equations  were  used  to  calculate  the 
reflectivity  spectra  for  the  VCSEL  structures  measured  and  designed  for  this  thesis  effort. 
The  program  used  was  developed  initially  at  the  University  of  Arizona.  The  algorithm  is 
designed  for  Al*Gai_*As  mirror  layers  and  cavity  materials,  but  should  work  reasonably 
well  for  cavities  containing  InGaAs,  which  has  a  refractive  index  reasonably  close  to  that 
of  GaAs.  The  algorithm  is  fairly  simple  and  is  explained  by  the  process  below: 

1.  Read  in  the  number  of  layers  in  the  top  and  bottom  mirrors  and  the  Aluminum 
concentrations  in  the  high  and  low  index  layers. 

2.  Read  in  the  length  of  the  cavity,  the  Aluminum  concentration  of  the  cavity,  and  the 
lasing  wavelength.  (Note:  The  Periodic  Gain  Structure  will  add  layers  to  the  internal 
structure  of  the  cavity.  For  this  analysis,  the  relatively  close  index  match  between 
InGaAs  and  GaAs  in  the  cavity  provides  a  negligible  effect  on  the  spectral  reflectance 
of  the  structure,  and  is  not  considered  in  the  analysis.) 

3.  Using  the  above  data  calculate  the  refractive  index  (rim)  via  the  above  equations. 
From  rim  aJid  tke  angle  of  incidence,  the  effective  index  {fim)  inay  be  determined 
using  equation  2.33. 

4.  Next,  calculate  the  phase  change  via  equation  2.34  above  for  each  successive  layer, 
then  determine  the  admittance  Y  (equation  2.40)  from  them;  Anally,  insert  Y  into 
equation  2.41  to  produce  the  desired  reflectivity. 


2-24 


Programs  to  calculate  the  reflectivity  of  a  single  multi-layer  stack  and  for  a  vertical 
cavity  structure  are  contained  in  Appendix  A. 

2.4-4  Spectral  Reflectivity  Stop  Band  MacLeod  recounts  the  derivation  for 
the  calculation  of  the  stop  band  for  the  quarter  wave  DBR  reflector  from  the  indices  of 
refraction  for  the  two  layers  in  the  stack.  The  relation  is  given  by  [26]: 


Au 

LJq 


4  .  .ifnH-riL 

—  sin  I - 

TT  KUh  +  IlL 


) 


(2.43) 


where: 

Aw:  stop  band  width  (frequency) 
wq:  center  frequency 

UHtTiL-  high  and  low  indices  for  the  layers 

From  equation  2.43,  it  is  evident  that  the  layer  thicknesses  do  not  affect  the  stop 
band  width,  and  that  changes  in  the  index  of  refraction  for  one  or  both  of  the  layers  can 
affect  it  considerably.  Variations  in  x  concentration  can  result  in  significant  changes  in 
the  indices  of  the  layers,  narrowing  or  broadening  the  stop  band.  As  a  result,  the  stop 
band  width  provides  an  effective  indication  of  the  control  of  x  concentration  in  the  layers, 
regardless  of  how  precise  the  control  on  layer  thickness. 

Equation  2.43  provides  a  rough  estimate  for  the  expected  stop  band  width  for  the 
structure.  In  practice,  dispersion  will  limit  its  effectiveness.  To  account  for  dispersive 
eflfects,  calculations  for  spectral  reflectivity,  and,  as  a  result,  stop  band  width,  should  be 
completed  via  the  matrix  method  described  earlier.  The  FORTRAN  programs  in  Appendix 
A  use  the  matrix  method. 


2.4.5  Penetration  Depth  The  penetration  of  the  optical  mode  into  the  mirrors 
will  affect  the  loss  (through  the  ttm/m  term)  in  equation  2.9,  and  the  effective  cavity  length. 
Both  of  these  concerns  are  addressed  in  this  section. 

Reflection  Delay  The  reflection  delay  Tr  is  used  to  express  the  change  in 
phase  of  the  optical  field  upon  reflection  from  the  mirror.  To  illustrate,  figure  2.8  shows 
an  equivalent  fixed  phase  mirror  model  of  the  cavity /mirror  interface  [4]. 
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VCSEL  Geometry  Equivalent  Phase  Geometry 

Figure  2.8.  Equivalent  Fixed  Phase  Mirror  Model  for  Determining  Reflection  Delay  in 
Multilayer  Stack  [After  Babic  [4]] 


For  the  linear  phase  mirror  at  the  top  of  figure  2.8,  the  observer  above  the  mirror 
surface  sees  a  reflected  optical  field  which  has  changed  in  amplitude  and  phase.  The 
amplitude  change  is  modeled  by  the  reflection  coeffUdent  (F^),  while  the  phase  change  is 
taken  into  account  by  the  (w  -  u;o)Tr  factor  in  the  exponential.  Here,  w  is  the  frequency 
of  the  incident  light,  while  (Jq  is  the  Bragg  frequency  (center  frequency  of  the  reflectivity 
stop  band).  The  resvdt  is  no  phase  change  for  center  Bragg  frequency  (u  =  Wq).  This  is 
consistent  with  the  mirror  design,  since  the  A/4  layers  are  designed  to  have  identical  phase 
(0  or  tt)  at  each  alternating  layer  at  the  Bragg  frequency.  Thus  the  light  entering  at  the 
Bragg  frequency  should  have  the  same  phase  when  it  exits  the  mirror.  Wavelengths  which 
are  removed  from  the  center  frequency  experience  the  linear  phase  shift  described  by  the 
frequency  difference  and  reflection  delay  in  the  exponential  ((w  —  u;o)rr). 

The  second  diagram  in  figure  2.8  models  the  linear  phase  mirror  as  a  fixed  phase 
mirror  set  back  from  the  observer’s  point  of  view  a  distance  Lr.  The  recovered  distance  (T,) 
is  the  penetration  depth  (sometimes  called  the  phase  penetration  depth).  The  refractive 
index  in  this  recovered  length  is  the  same  as  that  in  the  observer’s  medium. 

For  AJ*Gai_*As,  the  index  differences  in  the  high  (n^^)  and  low  (n^)  index  mirror 
layers  is  small  compared  to  the  average  index  of  the  layers  ( Off- n£  <<  na„).  The  reflection 
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delay  r,  for  small  index  differences  (uav  ss  n,  the  index  of  the  cavity)  has  been  obtained 
via  coupled  wave  theory  [23]  as: 


Tr 


n tanh  kI 

C  K 


(2.44) 


where: 

Tri  reflection  delay 
n:  index  of  refraction  of  the  cavity 
c:  speed  of  light  (free  space) 
k:  coupling  coefflcient 
/:  mirror  length 

The  coupling  coefficient  (/c)  for  square  wave  perturbations  (as  for  abrupt  index 
changes  in  mirror  layers)  is  k  =  2An/Ao.  Here  An  is  the  refractive  index  difference 
{nfj  —  ul)  in  the  mirror  layers,  and  Aq  is  the  Bragg  frequency. 


Phase  Penetration  Depth  Prom  figure  2.8  it  is  a  simple  exercise  to  de¬ 
termine  the  penetration  depth  and  fixed  phase  of  the  mirror.  The  phase  penetration  depth 
is  just  the  distance  the  light  travels  in  the  medium: 


. . .  and  the  fixed  phase  of  the  mirror  can  be  determined  by  equating  the  two  reflected 

fields. 


(2.45) 


...with  /?  =  nu/c  the  propagation  coefficient  of  the  mode.  Now  substitute  for  Lj 
and  solve  for  Gq: 


r  f,-Hu-Uo)Tr  _  p  p-HP2Lr-Bo) 
=>  (w  -  Ufo)Tr  =  2 - - - 00 

c  2n 
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UJQTr 


(2.46) 


=  ©0 


The  sum  of  the  phase  terms  in  the  equivalent  fixed  phase  mirror  cavity  can  be  equated 
to  an  integer  multiple  of  27r  to  determine  the  electromagnetic  modes  of  the  cavity  [4]: 


2P{Lfp  +  +  f'T/,,.)  —  ^oTrf^back)  “  ‘^0T’r(/rnl)  —  2jr/:  (2.47) 

where: 

Lpp:  design  length  of  the  cavity 
^T^ck.trni'  Phase  penetration  depth  in  the  back(front)  mirror 
Tr(haekjrnt)-  reflection  delay  for  the  back(front)  mirror 

...and  k  is  an  integer.  Note,  if  the  cavity  length  is  designed  to  the  same  center 
frequency  as  the  DBR  mirrors,  the  penetration  depth  and  reflection  delav  te^-ms  cancel,  so 
the  penetration  into  the  mirrors  does  not  affect  the  designed  resonant  mode  of  the  laser. 
Dispersion  will  affect  the  mode  spacing,  however.  The  penetration  will  be  affected  by 
dispersion,  since  the  change  in  refractive  indices  of  the  layers  alters  An.  This  will  alter  the 
longitudinal  mode  spacing. 

Mirror  Losses  Normally  the  mirror  losses  for  the  VCSEL  are  accounted 
for  by  the  penetration  depth  and  the  expected  loss  coefficient  for  the  material  (a^).  The 
energy  penetration  depth  is  usually  used  for  this  purpose.  Babic  and  Corzine  have  shown 
[4]  that  for  large  numbers  of  layers  (m  >  8)  the  phase  and  energy  penetration  depths 
approximate  each  other.  For  most  VCSEL  mirror  applications,  the  design  will  exceed  this 
constraint;  hence  one  can  use  the  phase  penetration  depth  as  the  mirror  penetration  depth: 


lm  =  Lr 


tanh  kI 
2k 


. .  .with  K  =  2An/Ao  as  defined  above. 


(2.48) 
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III.  Application  and  Performance 


3.1  Introduction 

The  application  of  the  equations  developed  in  chapter  2  to  the  building  of  a  vertical 
cavity  structure  is  the  aim  of  this  chapter.  Along  the  way,  the  equations  in  chapter  II 
will  allow  the  prediction  of  performance  in  the  VCSELs  to  be  grown.  The  initial  designs 
for  this  effort  focus  on  optically  pumped  VCSEL  structures.  Electrical  stimulation  of  the 
active  region  is  a  longer  range  goal  in  this  process.  A  summary  of  the  options  available  for 
electrically  stimulated  VCSELs  is  included  at  the  end  of  this  chapter. 

3.S  MQW  InGaAs  950  nm  Vertical  Cavity  Laser 

The  primary  constraints  on  the  cavity  design  for  this  thesis  were  the  lasing  wave¬ 
length  (950  nm)  and  the  mirror/cavity  materials.  The  materials  used  were  InrGai_eAs  for 
the  quantum  wells,  Al*Gai_,As  for  the  spacers,  barriers,  and  mirror  layers.  The  choice 
of  materials  was  dictated  by  the  available  materials  (and  number  of  cells)  for  the  MBE 
growth.  For  AlGaAs  quantum  well  VCSELs,  an  additional  aluminum  cell  would  have  been 
required  to  grow  two  different  Aluminum  concentrations  for  the  mirror  layers.  Rather 
than  reconfigure  to  grow  AlGaAs  quantum  well  VCSELs,  we  decided  to  pursue  InGaAs 
designs.  These  would  allow  the  mirrors  to  be  constructed  of  AlAs  (low  refractive  index 
layers)  and  GaAs  (high  refractive  index  layers).  The  chosen  VCSEL  wavelength  was  950 
nm,  but  initial  back  mirror  growths  were  designed  for  920  nm;  this  choice  was  dictated  by 
the  available  monochromator  (max  detectable  wavelength:  926  nm).  The  later  availability 
of  a  higher  range  monochromator  prompted  the  decision  to  pursue  a  950  nm  VCSEL,  along 
with  the  success  of  higher  wavelength  InGaAs  cavities  at  the  University  of  Arizona  [17] 
and  elsewhere  [20],  [16],  [8].  The  paragraphs  below  summarize  the  design  process  used  to 
construct  the  VCSELs  for  this  effort. 

3.2.1  Active  Medium  The  material  characteristics  for  In*Gai_,As  vary  with  x- 
concentration.  For  quantum  well  designs,  the  transition  energy,  and  as  a  result,  the  wave¬ 
length  of  the  emissions,  can  vary  with  the  band  gap,  the  effective  masses  of  the  holes  and 
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electrons,  and  the  well  width  (see  equation  2.31).  Many  combinations  of  these  parameters 
can  produce  the  desired  lasing  wavelength. 


S.2.2  InGaAs  Active  Region  Parameters  For  Inj;Gai_j,As  on  a  GaAs  sub¬ 
strate,  equations  for  the  energy  gap  and  the  effective  masses  of  the  electrons  and  holes  are 
[22]: 

m*  =>  me(i)  =  [.0665  —  0.0435i]  mo  (3-1) 

m„  ^  mhfc(x)  =  [0.45  —  0.04i]mo  (3-2) 

. .  .where  mo  is  the  free  electron  mass.  No  light-hole  effective  mass  equations  were 
available,  but  this  should  not  affect  the  estimate  greatly.  InGaAs  quantum  wells  on  GaAs 
substrates  undergo  only  compressive  strain  [12],  which  favors  the  heavy  hole  transition. 
The  band  gap  energy  for  bulk  In,Gai_,As  at  T=300K  is: 

£g(*)  =  1-424  -  1.614X  +  0.54x*  eV  (3.3) 

For  threshold  current  minimization,  the  material  gain  at  950  nm  as  a  function  of  x 
concentration  should  be  determined.  The  most  rigorous  procedure  would  be  to  plot  the 
spectral  gain  equation  [32]: 

where: 

g(E):  spectral  gain 
mo:  free  electron  mass 
n:  index  of  refraction 
Cq:  permittivity  of  free  space 
c:  speed  of  light  (vacuum) 

|M  |:  momentum  matrix  element  for  the  region 

Pr,d'  reduced  density  of  states 

fc  „:  Probability  band  state  is  occupied 
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The  use  of  this  equation  would  provide  the  material  gain  for  a  particular  frequency 
(E=  h  u)  which  could  be  plotted  vs  z-concentration  to  determine  the  concentration  that 
gives  the  maximum  material  gain.  To  do  so  it  is  necessary  to  know,  for  z-concentrations, 
the  electron  and  hole  effective  masses  (m',m„),  the  bulk  material  energy  gap  (Eg))  the 
spin  orbit  split-off  energy  (Aq),  and  the  energy  dispersion  of  the  conduction  band  (m*). 
These  parameters  have  been  determined  for  z-concentrations  for  AlrGai_*As  on  a  GaAs 
substrate  [2]  and  IniGai_eAsyPi_y  on  an  InP  substrate  [l],  [24],  but  not  all  have  been 
developed  for  Ini:Gai_,As  on  a  GaAs  substrate.  Even  without  these  parameters,  we  can 
still  make  an  informed  choice  for  the  Indium  (z)  concentration. 

Since  the  quantum  well  barriers  should  be  a  close  lattice  match  to  the  active  region, 
it  is  a  natural  choice  to  use  GaAs  for  the  barriers,  which  provides  a  close  lattice  match  to 
the  InGaAs  in  the  quantum  wells.  An  added  benefit  of  choosing  GaAs  for  the  barriers  is 
that  superlattices  with  GaAs  layers  have  evidenced  a  tendency  to  produce  very  good  layer 
uniformity  in  other  VCSEL  growth  reports  [14].  As  noted  in  chapter  2  the  barriers  should 
be  at  least  100  A  to  avoid  interwell  coupling.  With  the  considerations  above,  we  will  use 
100  A  thick  GaAs  barriers,  which  have  a  bulk  band  gap  of  1.424  eV  at  room  temperature 
[2]. 


3.2.3  Well  Width  and  z  Concentration  Now,  turning  to  the  transition  energy 
equation  (from  eq  2.31): 


/ii/  =  Eg(x)-H  + 

\m!(z) 


Selection  Rules  Yariv  notes  that  since  the  1=1  states  in  the  con¬ 
duction  and  valence  bands  have  the  highest  populations,  the  highest  optical  gain  is  to  be 
obtained  from  the  I  =  l{k  =  0)  conduction-to-valence  band  transition  (maximizes  (/«  -  /«) 
in  equation  3.4  above).  The  transition  energy  equation  then  reduces  to: 


hu  =  £„(x)  +  (^  + 


{/=l,fc  =  0}  (3.5) 
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A  plot  of  the  equation  above  for  A  =  950  nm  is  at  figure  3.1. 


Figure  3.1.  Required  Quantum  Well  Width  vs  Indium  Concentration  for  A  =  950  nm 
(1=1,  k=0) 


The  X  concentration  chosen  for  this  VCSEL  is  x  =  .2,  corresponding  to  a  64  Awell 
width.  In  choosing  this  well  width  it  is  necessary  to  consider  not  just  the  transition  energy: 

-  Barrier  Height:  This  width  provides  the  necessary  barrier  height:  the  bulk  en¬ 
ergy  gap  for  the  z  =  .2  concentration,  using  equation  3.3,  is  1.123  eV.  Corzine 
recommends  the  well-barrier  height  at  the  conduction  band  edge  be  at  least 
100  mev  [11]  to  preserve  the  energy  quantization.  Kolbas  [22]  says  the  conduc¬ 
tion/valence  band  split  is  70/30  for  InGaAs  on  GaAs.  The  gap  for  GaAs  (1.424 
eV)  should,  then,  be  large  enough  to  provide  the  necessary  well-barrier  height 
(.7  X  (1.424  -  1.123)  =  210mev). 

-  Critical  Thickness:  The  64  Awell  widths  should  be  safely  below  the  critical 
thickness  [15],  even  for  multiple  wells  [20].  The  barriers  will  be  100  Awide  to 
preclude  interwell  coupling,  and  this  is  also  below  the  critical  thickness. 

-  Fabrication  Capability:  The  64  Awell  is  within  the  capabilities  of  most  MB£ 
systems. 
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S.2.4  Cavity  Construction  Now  that  the  active  region’s  parameters  are  set,  the 
periodic  gain  structure  for  the  cavity  may  be  determined.  Cavity  length  is  determined  by 
the  lasing  wavelength  and  the  number  of  active  regions  to  be  included.  For  this  first-build, 
to  limit  growth  time  and  complications,  I  limited  the  number  of  quantum  wells  to  three. 
Additionally,  three  active  regions  were  chosen  for  the  PCS  structure.  Each  active  region 
having  three  quantum  wells  results  in  an  optical  thickness  per  region: 

tact  =  [3  X  6.4  X  Tlact  +  2  X  lOnjorr]  nTU 

With  Tiiarr  =  3.5455  the  index  of  refraction  for  GaAs  at  950  nm  using  equation  2.42 
outlined  on  page  2-23.  A  literature  search  turned  up  no  ready  means  of  determining  the 
index  of  refraction  for  varying  i-content  in  In^Gai.^As.  Fortunately,  the  refractive  index 
of  InAs  is  very  close  to  that  of  GaAs  (at  885  nm  nj„A,  =  3.696  and  UoaAs  =  3.633).  To 
estimate  the  index  for  Ino.2Cao.8As  we  interpolate  between  these  values  to  obtain  Uact  = 
3.5681.  These  indices  and  the  desired  well  and  barrier  thicknesses  correspond  to  an  optical 
thickness  tact  =  139.4  nm.  For  the  PGS  with  three  active  regions,  the  total  optical  thickness 
is  d  sa  418.25  nm.  The  cavity  was  to  be  designed  for  a  longitudinal  mode  of  950  nm.  If 
the  mirrors  are  designed  for  a  center  (Bragg)  frequency  at  950  nm,  we  can  eliminate  the 
penetration  depth  and  reflection  delay  terms,  as  discussed  in  the  Mirror  Design  section  of 
Chapter  II,  and  the  cavity  optical  length  can  be  determined  from  equation  2.47: 

^jiLpp  —  irrk 

Here,  13  =  27r/Ao,  Lpp  is  the  optical  length  of  the  cavity,  and  k  =  4,  since  we  want 
the  longitudinal  mode  to  peak  at  each  of  the  active  regions  (see  figure  3.2).  As  a  result, 
the  cavity  should  be  two  times  the  lasing  wavelength  or  a  /  =  1.9/zm  optical  thickness. 
This  corresponds  to  a  fill  factor  d/l  of  .22.  Corzine  recommends  fill  factors  in  the  vicinity 
of  0.2  -  0.3  for  optimal  confinement  factors  [11]. 

The  logical  choice  for  the  spacers  is  GaAs,  since  it  lattice  matches  well  to  both  the 
InGaAs  wells  and  the  AlGaAs  mirror  layers  adjacent  to  the  cavity,  and  provides  a  band 
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gap  sufficiently  removed  from  the  active  region  band  gap  to  minimize  absorption.  For  the 
parameters  discussed  so  far,  with  tact  =  139.4  nm  and  A/2  =  475  nm,  the  spacer  should  be 
335.58  nm  optical  thickness.  Since  the  spacer  is  GaAs  (n  =  3.546)  the  physical  thickness 
of  the  spacer  is  335.58/3.546  =  94.65  nm.  Now,  placing  the  active  regions  at  the  peaks  of 
the  electromagnetic  mode,  the  cavity  takes  the  form  of  figure  3.2.  The  spacers  adjacent 
to  the  mirrors  have  been  lengthened  so  that  the  three  active  regions,  two  interior  spacers, 
and  2  outside  spacers,  combine  for  a  total  cavity  optical  length  of  2  A. 

,  6tX 


Figure  3.2.  Cavity  for  the  Ino.jGao.sAs/GaAs  64  AQuantum  Well  VCSEL 


3.2.5  Mirror  Design  The  mirrors  for  the  cavity  design  from  the  previous  section 
were  designed  with  a  Bragg  frequency  of  950  nm.  The  use  of  InGaAs  for  the  active 
regions  allows  the  use  of  GaAs  in  the  mirrors,  since  the  higher  band  gap  of  the  GaAs 
mirror  layers  will  make  it  a  low  loss  material  at  the  lasing  wavelength.  The  mirrors  were 
constructed  of  Al,Gai_*As;  this  allowed  uninterrupted  growth  for  the  MBE  fabrication. 
The  natural  choice  for  the  concentrations  is  AlAs/GaAs.  This  combination  provides  the 
greatest  refractive  index  difference  for  the  AUGai_,As  family.  Large  An  is  important,  as 
the  number  of  layers  required  to  reach  a  desired  reflectivity  decreases  as  An  goes  up.  The 
VCSEL  will  be  grown  on  a  GaAs  substrate  to  provide  adequate  lattice  matching  for  the 
mirror  layers. 
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920  nm  Bragg  Frequency  Mirror  Design  The  initial  back  mirror  was 
grown  to  be  testable  with  the  reflectivity  measurement  systems  available.  At  920  nm  center 
frequency,  the  transmission  matrix  approach  outlined  in  Chapter  2  indicates  a  reflectivity 
of  99.17  %  can  be  achieved  for  the  mirror  with  17  periods  of  AlAs/GaAs  pairs,  starting 
with  AlAs  and  ending  with  GaAs.  The  parameters  for  the  mirror  are  below,  followed  by 
a  calculated  reflectivity  curve  in  figure  3.3,  produced  via  the  FORTRAN  program  DBR.f 
included  in  Appendix  A. 

Number  of  Layers:  34 

Substrate:  GaAs 

Peak  Reflectivity:  99.14  %  ®  920  nm  in  air 

Low  Index  Layers:  AlAs 
Calculated  Index:  2.98151 
Physical  Layer  width:  77.1  nm 

High  Index  Layers:  GaAs 
Calculated  Index:  3.5727 
Physical  Layer  width:  64.3  nm 
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Figure  3.3.  Bottom  Mirror  Reflectivity  as  seen  in  Air 

950  nm  Bragg  Frequency  Mirror  Design  The  back  mirror  for  the  final 
VCSEL  structure  will  be  grown  for  a  Bragg  frequency  of  950  nm.  At  950  nm  center 
frequency  the  transmission  matrix  approach  outlined  in  Chapter  2  indicates  a  reflectivity 
of  99.6  %  can  be  achieved  in  the  GaAs  spacer  with  19  1/2  periods  of  AlAs/GaAs  pairs. 
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starting  with  AlAs  and  ending  with  AlAs.  The  parameters  for  the  mirror  are  below, 
followed  by  a  calculated  reflectivity  curve  for  the  structure  as  seen  in  the  GaAs  spacer 
(figure  3.4),  produced  via  the  FORTRAN  program  DBR.f  included  in  Appendix  A. 

Number  of  Layers:  39 

Substrate:  GaAs 

Peah  Reflectivity:  99.6  %  @  950  nm  in  GaAs 

Low  Index  Layers:  AlAs 
Calculated  Index:  2.9742 
Physical  Layer  width:  79.85  nm 

High  Index  Layers:  GaAs 
Calculated  Index:  3.5455 
Physical  Layer  width:  66.99  nm 
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Figure  3.4.  Bottom  Mirror  Reflectivity  as  seen  in  GaAs  Spacer 

Top  Mirror  Design  The  VCSELs  to  be  designed  for  this  thesis  were  to 
have  the  output  beam  exiting  via  the  top  of  the  wafer.  With  this  in  consideration,  the 
bottom  mirror  design  is  straightforward;  it  requires  only  a  high  reflectivity  and  low  loss. 
The  top  mirror,  however,  must  provide  enough  reflectivity  to  satisfy  the  threshold  gain 
equation  (equation  2.3),  and  provide  a  low  enough  reflectivity  to  favor  the  beam  exiting 
the  top  mirror.  Most  successful  designs  have  resulted  after  several  iterations  of  fabrication 
and  measurements.  The  measurements  provide  an  accurate  means  of  determining  the 
gain/loss  profiles  of  the  materials  to  be  used  in  the  VCSEL  design.  For  a  first  cut,  Corzine 


[11]  has  achieved  success  with  a  geometric  mean  reflectivity  for  the  mirrors  of  >  95  %.  This 
reflectivity  is  on  the  order  of  those  used  in  successful  designs  at  the  University  of  Arizona 
[17],  and  at  AT&T  [19].  Given  the  bottom  mirror  reflectivity  (99.6  %)  and  this  minimum 
geometric  mean,  the  top  mirror  must  exceed  91  %  reflectivity.  To  provide  a  margin  of 
assurance,  I  chose  a  reflectivity  of  97  %.  This  top  mirror  reflectivity  provides  an  overall 
geometric  mean  reflectivity  of: 

R  =  \]RitmRtop  =  V.996  X  .97 
=  98.3% 


To  provide  the  top  mirror  reflectivity,  the  design  must  achieve  the  97  %  goal  with  the 
incident  medium  being  GaAs  and  the  exiting  medium  air.  Since  the  spacer  is  GaAs,  the 
first  layer  of  the  GaAs/AlAs  top  mirror  should  be  AlAs,  (a  GaAs  first  layer  next  to  the 
spacer  provides  no  added  reflection).  The  top  layer  should  be  GaAs,  since  it  provides  a 
larger  index  difference  with  the  air  interface,  and  GaAs  is  less  prone  to  oxidation.  As  a 
result,  we  will  have  an  even  number  of  layers  in  the  top  mirror.  To  determine  the  number 
of  layers,  we  use  equation  2.41  to  get  the  admittance  Y.  At  normal  incidence  fio  =  Wq; 
we  are  interested  in  the  reflectivity  from  the  spacer  looking  ou'.,  tiq  =  noaAt  ^  3.5455, 
”»u*  =  «oir  «  1.0  at  the  950  nm  center  wavelength: 


R  = 


Mo-y 
Mo  +  y 


^GoAt 

y  =  .0269 


-y 
+  Y 


With  this  admittance  Y,  and  knowing  we  need  an  even  number  of  layers  N,  we  can  use 
equation  2.40,  with  the  indices  of  refraction  for  GaAs  and  AlAs  at  normal  incidence  to 
determine  N. 


Y  = 


=> 


V  _  ^AlAt  ^0 

n^aA, 
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=»  .0269  = 


2.9742^1.0 

3.5455^ 


=;►  N  =  20S 
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Figure  3,5.  Top  Mirror  (22  layers)  Reflectivity  as  seen  from  GaAs  Cavity  Spacer 


To  exceed  the  97  %  goal,  we  choose  N  to  be  22.  The  reflectivity  curve  for  the  top 
mirror,  as  seen  from  the  GaAs  spacer  is  plotted  in  figure  3.5. 

The  overall  reflectivity  of  the  VCSEL  is  plotted  in  figure  3.6.  The  structure  has  a 
stop  band  of  125  nm  about  the  950  nm  center  frequency.  The  Fabry  Perot  resonzmce  is 
at  950  nm.  The  closest  minimum  is  at  883  nm,  this  is  the  likely  wavelength  location  for 
the  pumping  laser.  The  maximum  reflectivity  is  approximately  99.98  %;  air  being  the 
incident  medium.  Neglecting  dispersion  for  a  moment,  by  solving  equation  2.47  for  A:  =  3 
and  k  =  5  the  nearest  neighboring  modes  to  the  950  nm  resonant  frequency  would  be  1019 
nm  and  890  nm.  Both  modes  will  have  a  significantly  lower  gain  as  compared  to  the  950 
nm  longitudinal  mode  and  both  have  significantly  lower  reflectivity,  falling  at  the  edges  of 
the  stop  band.  These  factors  conspire  to  prevent  lasing  of  these  neighboring  longitudinal 
modes. 
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Figure  3.6.  Fabry- Perot  Effect  in  Designed  950  nm  VCSEL  using  39  layers  AlAs/HaAs 
bottom  mirror  and  22  GaAs/AlAs  layers  (top)  with  2  A  MQW  spacer 

3.3  Carrier  Guided  VCSELs 

The  VCSELs  measured  in  this  research  were  not  processed  further  after  their  MBE 
growth.  It  is  common  to  impart  a  mechanism  to  provide  a  means  of  guiding  the  carriers 
which  will  be  optically,  or  electrically,  injected.  Gain  guiding  and  index  guiding  are  the 
typical  choices. 


Coin  CuuUof  Index  Guiding 

Figure  3.7.  Gain  Guiding  and  Index  Guiding  in  Conventional  Processed  VCSEL  Struc¬ 
tures 

Designed  gain  guiding  is  usually  accomplished  by  ion  implantation  on  the  wafer  (see 
figure  3.7).  The  ion  implantation  is  achieved  by  conventional  masking  and  implantation 
techniques  [19].  Regions  are  formed  in  the  wafer  which  inhibit  gain  in  volumes  adjacent  to 
the  volume  where  lasing  is  desired.  The  relatively  high  gain  region  tends  to  concentrate 
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the  carrier  population  in  the  desired  volume,  favoring  stimulated  emission.  Index  guiding 
is  less  subtle.  By  etching  the  mirror  and  cavity  layers  away  from  the  substrate  to  form  a 
columnar  structure,  the  air/semiconductor  index  difference  at  the  sides  of  the  column  is 
used  to  contain  the  carriers. 

The  VCSELs  tested  for  this  effort  did  not  use  ion  implantation  or  etching  to  provide 
a  guide  mechanism  for  the  carriers.  With  sufficient  pump  power  and  a  small  enough  pump 
spot,  the  optically  generated  carriers  can  form  their  own  guide  mechanisms.  The  injection 
of  the  carriers  tends  to  elevate  the  gain  in  the  area  around  the  pump  spot;  in  addition, 
the  high  concentration  of  carriers  also  tends  to  increase  the  local  refractive  index  of  the 
medium.  Thus,  the  carriers  form  a  self  generated  gain  and  index  guiding  mechanism  to 
favor  the  population  inversion  [13].  This  type  of  carrier  induced  guiding  is  illustrated  in 
figure  3.8. 


Active  Region 


Carrier  Generated  Gain/Index  Guiding 
Figure  3.8.  Gain  and  Index  Guiding  from  Carrier  Injection 


3.4  Iigection  Lasers:  Electrical  Pumping 

The  design  so  far  has  focused  on  the  building  of  optically  pumped  surface  emitters. 
The  formation  of  a  junction  in  the  VCSEL  doesn’t  present  much  difficulty,  the  entire 
structure  is  of  semiconductor  and  it  is  relatively  straightforward  to  impart  dopants  to  the 
layers  as  necessary  to  form  the  desired  acceptor/donor  levels.  A  p-i-n  structure  is  the 
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commonly  used  scheme,  leaving  the  active  region  undoped  to  lower  cavity  losses  [19].  The 
difficulty  in  electrical  injection  is  that  the  same  abrupt  discontinuities  in  the  mirror  layers 
that  provide  efficient  DBR  reflection  result  in  high  resistivity  of  the  structure.  As  with 
other  advances  in  VCSEL  design,  the  lead  of  in-plane  semiconductor  lasers  provided  a 
guide  for  early  VCSEL  designs  in  this  area. 

S.4‘1  Active  Region  Design  A  natural  improvement  to  the  abrupt  well-barrier 
structure  of  the  quantum  wells  used  for  this  thesis  is  the  Graded  Index  Separate  Con¬ 
finement  Heterostructure  (GRIN  SCH).  The  GRIN  SCH  employs  a  superlattice  at  the 
well-barrier  interface  to  provide  a  smoother  transition  from  the  well  to  the  barrier.  The 
separate  confinement  heterostructure  increases  the  confinement  factor  (F),  which  falls  off 
as  the  square  of  the  thickness  of  the  active  region  (~  t^)  [32].  As  the  quantum  wells  be¬ 
come  thinner,  the  reduction  in  F  can  become  significant.  The  GRIN  SCH  is  more  effective 
than  abrupt  interfaces  at  capturing  carriers  in  the  wells.  This  becomes  important  for  thin 
quantum  wells. 

There  is  a  residual  benefit  of  this  structure  for  current  injection  purposes.  The 
gradual  variations  in  the  energy  gap  (compared  to  those  for  abrupt  interfaces)  between 
materials  provide  less  dramatic  well/barrier  discontinuities.  This  lowers  resistance  of  the 
structure  [19].  Non-quantum  well  structures  can  also  use  graded  index  interfaces  to  reduce 
resistances  in  the  active  regions.  The  application  is  a  straightforward  extension  of  that 
used  in  the  injection  design  for  the  DBR  mirrors  described  below.  The  number  of  interfaces 
to  be  concerned  about  in  either  UGS  or  PGS  active  regions  is  relatively  few  compared  to 
the  DBR  mirrors.  As  a  result,  the  pay  back  for  the  design  effort  is  less  dramatic. 

3.4-S  DBR  Mirror  Design  The  abrupt  discontinuities  of  the  DBR  mirror  as 
described  for  the  optically  pumped  structures  of  this  thesis  will  provide  relatively  high 
resistivities  if  used  in  an  injection  scheme.  In  the  p-type  layers,  the  periodic  well-barrier 
discontinuities  produce  the  high  resistances.  In  the  n-type  layers,  the  high  resistance  is 
attributed  to  poor  conductivity  of  the  AlAs  layers  [19]. 
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The  remedies  to  the  resistance  of  the  mirror  layers  include  softening  the  energy  gap 
transitions  (non-abrupt  interfaces)  or  etching  through  the  mirrors  to  inject  carriers  directly, 
avoiding  the  discontinuities  in  the  mirror  layers.  These  two  methods  are  shown  in  figure 
3.9  below. 


Noo-Abnipt  Interfaces  (Supcrlaitice)  Tnmsvase  Cunent  Flow 

Figure  3.9.  Two  Methods  of  Avoiding  High  Mirror  Layer  Resistances  for  Current  Injec¬ 
tion  in  VCSELs 


Non-Abrupt  Interfaces  The  use  of  non-abrupt  interfaces  has  met  very 
good  success.  Two  methods  are  prevalent.  Th^  first  is  to  limit  the  x  concentrations  of  the 
two  mirror  layers  so  that  their  difference  is  less  than  Ax  =  0.6.  This  approach,  however, 
increases  light  penetration  into  the  mirrors,  which  will  increase  mirror  loss,  and  require 
more  layers  to  achieve  a  desired  reflectivity. 

The  second  method  (which  is  shown  in  figure  3.9)  is  to  incorporate  superlattice 
transitions  between  the  layers  in  the  DBR  mirrors.  Tai  achieved  two  orders-of-magnitude 
reduction  in  resistivity  using  a  10  A,  10  period  superlattice  at  each  interface.  The  super 
lattice  was  Alo.7Gao.3As  /GaAs  and  the  quarter  wave  layers  were  also  Alo.7Gao.3As  /GaAs 
[31]. 
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Transverse  Current  Flow  As  a  promising  alternative  to  non-abrupt  in¬ 
terfaces,  Jewell  suggests  the  use  of  a  transverse  current  in  high  conducting  layers  above 
and  below  the  active  region.  The  current  would  be  injected  by  etching  through  the  mirrors 
to  place  contacts.  This  method  stems  from  the  power  efficiency  equation  of  the  VCSEL 
structure,  given  by  [19]: 


Pejf  —  Veff  X 


V,  -1- 


(3.6) 


where: 

Pejj-  power  efficiency 
rjejj:  differential  quantum  efficiency 
/op:  operating  current 
Itk'-  threshold  current 
Vf’.  normalized  bandgap  energy 
R:  series  resistance  of  the  VCSEL  structure 

This  equation  shows  how  an  increased  voltage  drop  (lopR)  lowers  the  power  efficiency. 
To  lower  series  resistance  R,  vertical,  rather  than  transverse,  current  flow  is  desired  through 
the  active  region.  Transverse  current  flow  through  the  undoped  active  region  increases  the 
voltage  drop,  which  decreases  efficiency  High  conductivity  layers  would  be  used  above  and 
below  the  active  region  for  transverse  current  flow,  and  current  flow  through  the  active 
region  would  then  be  vertical.  The  drawback  of  this  structure  is  the  extra  processing 
required  to  produce  the  desired  current  flow.  This  extra  processing  is  compounded  as  the 
packing  density  is  increased. 
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IV.  Experimental  Method 


The  experimental  part  of  this  research  involved  measuring  output  beam  character¬ 
istics  and  surface  reflectivities.  All  of  the  output  beam  characteristics  and  many  of  the 
surface  reflectivity  measurements  were  performed  using  one  experiment  configuration.  Ad¬ 
ditional  reflectivity  measurements  with  an  alternate  configuration  were  performed  to  com¬ 
pensate  for  some  of  the  equipment  limitations  encountered  in  covering  the  wavelength 
range  for  the  InGaAs  structures. 

4-1  Laboratory  Configuration 

The  equipment  configuration  used  to  optically  pump  the  samples  is  shown  in  figure 
4.1.  In  this  configuration,  the  pump  beam  enters  the  VCSEL  vertically  via  the  top  mirror. 
A  short  focal  length  lens  is  used  to  concentrate  the  pump  onto  a  small  area  of  the  structure 
(see  figure  4.2). 


Figure  4.1.  Equipment  Configuration  used  for  Output  Beam  Characterization 


The  output  beam  exits  via  the  top  of  the  wafer  for  all  of  the  samples  measured  in  this 
experiment.  The  output  beam  travels  back  down  the  same  optical  path  as  the  reflected 
pump  beam.  The  pump  and  output  are  then  separated  via  filters  and  beamsplitters  and 
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directed  to  the  measurement  devices.  A  brief  description  of  the  major  elements  in  the 
configuration  foUows,  Table  4.1  provides  relevant  specifications  for  the  equipment: 


VCSEL 


I .  j 

Surface  ^  ^ 


Continuously 
variable  oblique 
incidence 


Figure  4.2.  Optical  Pump  Beam  and  Focusing  Lens  Geometry 
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Ar"*"  Laser 


Dye  Laser 


Dichroic  Beamsplitter 


Focusing  Lens 

Positioner 

Power  Meter 

Filter 

TEC 


Monochromator 


Provides  the  means  to  optica.Jy  pump  the  dye  laser.  The  Argon 
laser  is  the  optical  power  source  for  the  experiment 

Chosen  to  provide  a  tunable  pumping  source.  The  dye  used  for 
this  application  (LDS-821)  provided  a  wavelength  range  (790 
nm  -  905  nm)  which  covered  the  typical  pump  wavelengths  for 
both  GaAs  and  inGaAs  VCSELs.  The  tunable  dye  laser  also 
provided  a  means  of  measuring  the  spectral  reflectivity  of  the 
samples. 

The  DBS  used  to  pump  the  GaAs  VCSELs  had  a  maximum 
transmissivity  at  830  nm  and  a  maximum  reflectivity  at  870  nm. 
This  property  allowed  the  incident  pump  to  pass  from  the  dye 
laser.  After  reaching  the  sample,  the  output  laser  and  the  pump 
return  to  the  DBS  (see  figure  4.1).  The  DBS  transmits  most  of 
the  pump  and  reflects  most  of  the  output  beam.  This  provides 
.part  of  the  filtering  of  the  pump  from  the  VCSEL  output. 

Short  focal  length  (8  mm)  lens  which  focused  the  pump  beam 
on  the  wafer. 

An  accurate  positioner  was  necessary  to  provide  repeatable  po¬ 
sitioning  of  the  pump  onto  the  wafer  surface. 

used  to  measure  pump  power  in,  reflected  power  from  the  wafer, 
and  VCSEL  power  out. 

was  used  to  filter  the  remainder  of  the  reflected  pump  from  the 
VCSEL  output  before  allowing  the  VCSEL  component  to  pass 
to  the  measurement  devices. 

The  Thermoelectric  cooler  (TEC)  was  used  to  maintain  constant 
temperature  for  the  samples,  which  tended  to  heat  from  the 
incident  pump  energy.  The  TEC  was  mounted  between  copper 
plates.  The  bottom  copper  plate  was  attached  to  a  heat  sink 
which  was  later  glued  to  the  top  of  the  positioner.  The  VCSEL 
was  bonded  with  heat  sink  compound  to  the  top  copper  plate. 
All  surfaces  between  the  VCSEL  wafer  and  the  heat  sink  were 
coated  with  heat  sink  compound  to  form  thermal  bonds. 

The  monochromator  was  used  to  determine  the  output  wave¬ 
length  of  the  VCSEL  as  well  as  the  pump  wavelength  (when  the 
narrowband  filter  is  removed).  A  detector  was  used  to  monitor 
the  exit  slit  of  the  monochromator. 
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IR  Camera  The  pump  and  VCSEL  output  were  separately  directed  to  the  IR 
camera  to  image  the  spots  from  the  top  of  the  wafer.  This  allowed  a 
measure  of  the  laser  spot  diameter,  by  comparison  to  the  size  of  the 
imaged  pump  spot.  The  minimum  spot  size  for  the  focusing  lens  at 
the  VCSEL  provided  the  pump  spot  size.  The  image  of  the  pump 
spot  also  allowed  a  qualitative  assessment  of  the  wafer  surface. 

Polarizer  Cube  Used  to  determine  the  polarization  of  the  output  beam.  Component 
polarizations  were  determined  by  inserting  a  quarter  wave  plate  to 
rotate  the  light  incident  from  the  sample. 


Table  4.1.  Equipment  Specifications  for  Laboratory  Configuration 


Ar"*"  Laser 

Model 

Spectra-Physics  2020-03 

Wavelength 

514.5  nm 

Output  Power 

Dye  Laser 

Model 

Spectra-Physics  375-50 

Dye 

LDS-821  (Styril  9) 

Wavelength 

Tunable  790  -  905  nm 

Output  Power 

~  200  mW  (typical) 

Linewidth 

~  666  A 

Beam  Divergence 

1.6  mrad 

Beam  Mode 

TEMoo  Vertically  Polarized 

Mode  Spacing 

420  MHz 

Dichroic  Beam 

Mfr 

CVI 

Splitter 

lUno* 

870  nm 

I'maT 

830  nm  at  0  =  5® 

Filter 

Mfr 

Ealing 

Center  Wavelength 

880  nm 

BW  (FWHM) 

10  nm 

Monochromator 

Model 

Ebert  82-546 

Wavelength  Range 

200  nm  -  926  nm 

Accuracy 

±  2  A 

Power  Meter 

Model 

Coherent  Fieldmaster 

Detector 

LM-2  Silicon  Sensor 

Focusing  Lens 

Mfr 

Melles  Griot 

fo 

8  mm 

spot  size 

1.02  fim 

working  distance 

1.13  mm 
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4-1.1  Dye  Laser  Configuration  Refiectivity  Measurements  Reflectivity  mea¬ 
surements  using  the  configuration  in  figure  4.1  were  performed  with  and  without  the  fo¬ 
cusing  lens.  The  lens  provided  a  means  of  obtaining  the  reflectivity  for  a  small  spot  on 
the  wafer.  This  allowed  reflectivity  measurements  for  adjacent  lasing  and  non-lasing  spots 
to  explore  the  link  between  reflectivity  and  lasing  behavior.  The  reflectivities  obtained 
with  the  lens  in  place  are  difficult  to  predict  via  the  theory  summarized  in  chapter  II.  The 
calculation  is  complicated  by  the  action  of  the  lens  which  focuses  the  collimated  beam 
onto  thf  spot  (see  figure  4.2).  The  light  impinging  on  the  spot  is  no  longer  normal  to  the 
surface.  The  pump  rays  are  at  varying  degrees  of  oblique  incidence,  depending  on  where 
they  enter  the  lens.  The  resultant  reflectivity  is  a  convolution  of  the  obliquely  incident  ray 
reflectivities. 

The  lens  measurements  are  more  productive  as  a  comparison  from  spot  to  spot  on  the 
wafer,  rather  than  as  a  comparison  to  theoretical  calculations.  Additional  measurements 
were  done  on  the  wafers  without  the  focusing  lens  to  compare  measured  to  predicted 
reflectivity  for  orthogonal  incidence.  Measurements  without  the  lens  can  provide  a  means 
of  determining  the  quality  of  the  fabrication  of  the  wafer.  For  example,  systematic  thickness 
variations  show  as  a  shift  in  the  reflectivity  spectrum  by  changing  in  equation  2.34.  A 
ch2mge  in  the  Aluminum  concentration  of  one  of  the  layers  will  change  the  width  of  the 
stop  band,  as  noted  in  section  2.4.4.  The  drawback  to  the  measurements  taken  without 
the  focusing  lens  is  that  the  spoc  on  the  wafer  is  much  larger  (  ~  3mm  ),  so  the  reflectivity 
for  the  area  is  averaged  over  the  <irea  of  the  spot,  which  takes  in  both  lasing  and  nonlasing 
locations. 

4-1-2  Additional  Refiectivity  Measurements  These  measurements  were  per¬ 
formed  to  measure  a  broader  spectral  range  than  could  be  obtained  from  the  measurement 
configuration  in  figure  4.1.  The  configuration  is  shown  in  figure  4.3,  and  consists  of  an 
Argon  ion  laser  pumping  a  tunable  Ti-Sapphire  laser.  The  same  power  meter  used  for 
the  dye  laser  measurements  was  used  for  these  added  reliectivity  measurements.  The  Ti- 
Sapphire  configuration  provided  a  much  broader  measurement  range  (850  nm  -  1.0  fim  ) 
for  the  reflectivity  measurements. 
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These  measurements  were  performed  without  a  focusing  lens  for  the  incident  beam, 
and  represent  the  reflectivity  for  a  broad  (  ~  3mm  )  spot  on  the  wafer.  To  meausure  the 
reflected  beam  off  the  sample,  the  incident  beam  could  not  be  orthogonal  to  the  sample 
Oblique  incidence  shifts  the  reflectivity  spectrally,  as  discussed  in  the  Mirror  section  of 
chapter  II.  The  measured  reflectivities  provided  in  chapter  V  are  plotted  against  oblique 
incidence  theoretical  curves  to  compensate. 
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V.  Experimental  Results 


The  configuration  described  in  chapter  IV  was  used  to  measure  characteristics  from 
two  VCSEL  wafers  with  GaAs  active  regions.  These  samples  were  grown  at  the  University 
of  Arizona.  Additionally,  measurements  on  the  back  mirror  for  the  VCSEL  designed  for 
this  thesis  were  completed.  Scheduling  difficulties  prevented  completion  of  the  full  VCSEL 
structure  designed  in  chapter  III.  The  theoretical  curves  presented  here  are  based  on  the 
theory  outlined  in  chapters  II  and  III,  and  used  the  programs  in  Appendix  A  to  perform 
the  calculations. 

5.1  Measurement  Parameters 

The  parameters  measured  for  each  sample  were  lasing  wavelength,  reflectivity,  po¬ 
larization,  lasing  mode,  and  pump  wavelength. 

5.1.1  Reflectivity  Each  sample  has  a  predicted  reflectivity  curve  based  on  the  growth 
parameters  of  the  sample.  These  curves  were  calculated  using  the  FORTRAN  programs 
in  Appendix  A.  Both  lasing  and  non-lasing  locations  on  the  samples  were  measured  for 
reflectivity  vs  wavelength.  The  measured  wavelength  range  bracketed  the  upslope  of  the 
Fabry-Perot  effect  for  the  sample.  Loss  in  the  mirror  and  cavity  layers  due  to  absorption 
and  scattering  was  not  factored  in  for  these  curves.  As  a  result,  the  theoretical  spectral 
reflectivity  curves  can  be  much  higher  in  magnitude  than  the  measured  curves. 

The  pump  laser  (dye  laser)  was  attenuated  to  a  level  below  that  which  would  start 
VCSEL  lasing;  so  that  only  surface  reflectivity  was  collected.  Input  and  output  power 
were  measured  at  the  points  indicated  in  figure  4.1.  The  reference  reflectivity  for  the 
samples  was  a  high  reflectivity  broadband  mirror,  substituted  for  the  VCSEL  sample.  The 
designed  reflectivity  of  this  mirror  was  >  99.5  %  over  the  wavelength  range  measured.  The 
optical  path  for  both  the  reference  and  the  VCSEL  sample  included  the  turning  mirror, 
focusing  lens,  and  dichroic  beamsplitter  losses  (see  figures  4.1  and  4.2).  The  input  to  output 
power  ratio  for  the  VCSEL  sample  was  then  compared  to  this  reference  to  determine  the 
reflectivity  of  the  sample. 
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5.1.2  Lasing  Wavelength  Lasing  wavelength  was  measured  using  the  monochro¬ 
mator  indicated  in  Table  4.1.  The  wavelength  accuracy  of  the  monochromator  was  ±  2 

A. 


5.1.3  Pump  Wavelength  Indicated  the  location  of  the  pump  that  favored  highest 
lasing  output  power.  This  characteristic  was  measured  with  the  same  monochromator  used 
to  measure  the  lasing  wavelength. 

5.1.4  Polarization  The  VCSELs  measured  in  this  research  had  GaAs  bulk  active 
regions  (UGS).  The  close  lattice  match  of  these  active  regions  to  the  Al;fGai_xAs  mirror 
layers  provides  very  little  strain.  As  a  result,  random  polarization  was  expected.  The 
VCSEL  samples  used  were  small  pieces  of  the  original  wafer  with  no  growth  directions 
indicated.  Without  a  growth  direction,  the  choice  of  which  output  polarization  is  horizon¬ 
tal,  or  vertical,  is  arbitrary.  I  chose  to  relate  the  polarization  of  the  output  of  the  VCSEL 
to  the  input  pump  from  the  dye  laser.  The  dye  laser  output  is  vertically  polarized  (see 
table  4.1),The  dye  laser  input  polarization  should  not  affect  the  VCSEL  output  polariza¬ 
tion,  since  once  the  carriers  are  injected,  stimulated  emission  in  the  GaAs  active  region 
forms  its  own  polarization  direction.  The  reflected  dye  laser  pump,  then,  serves  only  as  a 
convenient  measurement  reference  for  the  VCSEL  output  polarization. 

All  polarizations  measured  were  linear  and  in  the  same  direction.  Rotation  of  the 
dye  laser  pump  beam  polarization  did  not  affect  VCSEL  output  polarization.  Since  only 
four  locations  on  each  wafer  were  measured,  it  would  be  premature  to  draw  conclusions 
concerning  the  polarization  tendencies  of  these  structures.  The  capability  to  measure 
the  VCSEL  output  polarization,  however,  is  now  available  for  future  studies  with  this 
configuration. 

5.1.5  Mode  A  check  on  the  lasing  spot  imaged  in  the  infrared  camera  provided  this 
measurement.  When  higher  order  transverse  modes  were  encountered,  input  power  and 
pump  beam  positioning  were  varied  to  see  if  they  could  alter  the  modes. 
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5.2  Sample  1:  3  X  UGS  GaAs  Active  Region 

Sample  1  consisted  of  a  3  A  long  gain  medium  of  GaAs,  the  3  X  sample  was  the  oldest 
sample  used  for  this  study;  it  was  grown  at  the  University  of  Arizona  in  November  of  1990. 
The  sample  geometry  is  represented  in  figure  5.1  below,  along  with  the  locations  on  the 
surface  where  measurements  were  taken  for  lasing  and  non-lasing  spots. 


Figure  5.1.  Sample  1:  3  A  Active  Region  VCSEL  Geometry 


5.2.1  Material  Parameters  The  material  parameters  for  this  sample  are  summa¬ 
rized  below: 


Number  of  Layers  (front): 

34 

Number  of  Layers  (back): 

45 

Active  Region: 

GaAs  UGS 

Active  region  thickness 

.726  ;im 

Substrate: 

GaAs 

Design  wavelength: 

875  nm 

Low  Index  Layers: 

AlAs 

Calculated  Index  @  875  nm: 

2.99573 

Physical  Layer  width: 

73.2  nm 

High  Index  Layers: 

Al.i2rGa,873. 

Calculated  Index  ^  875  nm: 

3.52773 

Physical  Layer  width: 

61.8  nm 

The  theoretical  reflectivity  of  Sample  1  based  on  the  parameters  described  above 
(with  no  loss  assumed)  is  represented  in  figure  5.2.  There  is  a  slight  dip  in  the  curve  at 
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the  design  wavelength  for  the  structure,  due  to  the  Fabry- Perot  effect  of  the  mirror/spacer 
integrated  etalon.  The  minimum  nearest  the  stop  band  is  at  w  SSOnm 
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Figure  5.2.  Sample  1:  Theoretical  Spectral  Reflectivity  near  875  nm 


5.2.2  Sample  1  Measurements  As  shown  in  figure  5.1,  lasing  spots  for  sample 
1  were  grouped  in  the  lower  left  quadrant  of  the  sample.  Reflectivity  measurements  for 
two  of  the  spots,  and  output  beam  measurements  for  all  of  the  lasing  spots  follow.  The 
reflectivity  measurements  for  two  lasing  spots  are  paired  with  reflectivity  measurements 
for  neighboring  non-lasing  spots.  Additional  reflectivity  measurements  at  non  lasing  areas 
on  the  sample  are  also  included. 

Sample  1;  Spot  1  : 

Lasing  Wavelength  872.0  nm 
Pump  Wavelength  833.0  nm 
Polarization:  Horizontal 
Mode:  TEMqo 

Reflectivity  The  reflectivities  measured  for  spot  1  indicate  that  the 
lasing  and  neighboring  non-lasing  areas  on  the  wafer  surface  have  very  nearly  the  same 
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reflectivities.  This  indicates  that  the  Aluminum  concentrations  and  layer  thicknesses  are 
fairly  uniform.  These  reflectivities  represent  a  convolution  of  oblique  incidence  beams,  the 
result  of  the  use  of  the  focusing  lens  (see  figure  4.2)  to  pump  the  VCSEL.  The  location  of 
the  pump  minimum  at  833  nm  indicates  the  structure  has  a  reflectivity  near  the  design 
minimum.  Oblique  incidence  shifts  the  spectral  reflectivity  to  shorter  wavelengths;  the 
lens  convolves  the  incident  beams  to  average  the  total  shift.  This  indicates  the  structure 
may  have  slightly  thicker  mirror  or  spacer  layers  than  designed. 

The  measured  reflectivities  do  not  reach  the  99  %  mark  expected  for  the  lower  wave¬ 
lengths  of  the  stop  band.  Absorption  and  scattering  losses,  which  were  neglected  in  the 
calculated  spectral  reflectivity  curve,  contribute  to  the  lower  reflectivity. 
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Figure  5.3.  Sample  1  Measured  Spectral  Reflectivities  for  Lasing  and  Non-Lasing  Loca¬ 
tions  near  Spot  1 


Sample  1;  Spot  2  : 

Lasing  Wavelength  872.1  nm 
Pump  Wavelength  830.9  nm 
Polarization:  Horizontal 
Mode:  TEMoo 
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Reflectivity  The  reflectivities  measured  for  spot  2  are  also  very 
similar  from  the  lasing  to  non-lasing  areas.  The  lasing  spot’s  reflectivity  curve  in  this  case 
is  distinctly  higher  at  the  longer  wavelengths.  The  first  minimum,  however  is  nearly  the 
same  magnitude  as  for  the  non-lasing  spot. 

The  difference  in  reflectivities  for  the  wavelengths  in  the  stop  band  are  greater  for 
spot  2.  In  all  of  the  lasing/non-lasing  areas,  the  imaged  surface  of  the  VCSEL  wafer  at 
830  nm  indicated  a  small  (~  l/xm  as  compared  to  the  pump  spot)  ‘island’  on  the  wafer. 
When  the  pump  spot  was  focused  on  this  ‘island’,  lasing  would  occur;  movement  of  the 
pump  off  the  ‘island’  ceased  lasing.  The  effect  could  not  be  captured  (photographed)  with 
the  configuration  used  for  this  experiment.  The  ability  to  see  (in  the  imaged  spot)  the 
physical  difference  indicated  there  might  be  a  change  in  reflectivities  of  the  two  areas.  In 
some  (such  as  spot  2  here)  the  difference  in  reflectivity  was  more  dramatic. 
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Figure  5.4.  Sample  1  Measured  Spectral  Reflectivities  for  Lasing  and  Non-Lasing  Loca¬ 
tions  near  Spot  2 


Sample  1;  Spot  3  and  Spot  4  Only  lasing  data  are  presented  for  Spot  3 
and  4,  as  the  reflectivity  data  for  these  spots  is  similar  to  that  for  Spots  1  and  2,  and  is 
not  necessary  for  the  subsequent  analysis. 
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Spot  3: 


Lasing  Wavelength  874.1  nm 
Pump  Wavelength  832.0  nm 
Polarization:  Horizontal 
Mode:  TEMoo 

Spot  4: 

Lasing  Wavelength  880.7  nm 
Pump  Wavelength  834.0  nm 
Polarization:  Horizontal 
Mode:  TEMoo 

Sample  1;  I^on-Lasing  Spots  :  Spots  5,  6  and  7  in  figure  5.1  did  not 
lase.  In  fact,  they  were  representative  of  the  areas  on  the  wafer  which  did  not  even  start 
photoluminescence.  The  reflectivity  data  for  these  areas  is  contained  below  and  is  included 
to  allow  comparison  to  the  lasing  areas  for  Sample  1  already  presented. 


Sample  1 ,  Reflectivities  for  Non-Lasing  Spou  3,6,  and  7 
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Figure  5.5,  Sample  1  Measured  Spectral  Reflectivity  for  Non-Lasing  Locations:  Spots  5, 
6,  and  7 
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5.2.3  Sample  1  Analysis  The  lasing  wavelengths  for  this  sample  indicate  the 
lasing  behavior  matched  the  design  expectations  fairly  closely  (A  «  875nm).  Spots  1,  3, 
and  4  indicated  pump  locations  near  the  first  dip  in  the  Fabry-Perot  curve  adjacent  to 
the  stop  band.  A  slight  shift  to  shorter  wavelengths  could  be  due  to  the  oblique  incidence 
forced  by  the  focusing  lens. 

Polarization  was  linear  for  all  measured  lasing  spots  on  the  wafer.  The  polarization 
direction  was  the  same  for  all  spots.  Jewell  indicates  in  a  1989  study  that  polarization  did 
not  favor  a  particular  direction  for  GaAs  VCSELs  [21].  The  spots  measured  here  indicate 
a  possible  tendency;  however,  larger  samples,  with  more  lasing  spots,  and  more  complete 
details  about  the  growth  directions  are  needed  to  draw  conclusions  concerning  polarization. 

The  mode  for  all  four  spots  on  sample  1  was  TEMoo-  A  picture  captured  from  spot 
1  is  shown  in  figure  5.6.  The  pump  spot  shown  is  imaged.  The  pump  spot  was  assumed 
to  be  the  spot  size  for  the  diode  laser  lens  used  to  focus  the  collimated  dye  laser  beam. 
The  specified  spot  size  for  collimated  input  for  this  lens  is  1.02  /xm  (see  table  4.1).  This 
provides  a  measure  of  the  spot  size  for  the  VCSEL  output  relative  to  the  «  1.02/im  pump 
spot;  the  VCSEL  spot,  then  is  approximately  2.5  /xm  in  diameter. 


Figure  5.6.  Sample  1,  Spot  1  :  TEMoo  Lasing  and  Pump  Spots 


The  reflectivity  curves  for  each  spot  are  nearly  identical  (see  figure  5.3  and  5.4). 
With  this  close  match,  it  is  unlikely  that  spectral  reflectivity  alone  can  be  used  to  predict 
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lasing  behavior  in  VCSELs.  Spots  5,  6,  and  7,  which  were  from  an  area  of  the  VCSEL 
wafer  which  evidenced  no  tendency  to  lase,  show  low  reflectivities  where  the  stop  band 
should  be  («  80%),  indicating  reflectivity  measurements  could  be  used  as  a  cursory  check 
on  the  wafer’s  tendency  toward  lasing.  However,  if  pursuing  a  quick  means  of  determining 
fabrication  success,  final  conclusions  about  lasing  quality  and  the  abundance  of  lasing  spots 
on  the  sample  would  require  a  more  than  reflectivity  measurements. 
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5.3  Sample  2:  3  A  UGS  GaAs  Active  Region 

Sample  2  had  a  design  wavelength  similar  to  that  of  sample  1,  but  had  a  greater 
number  of  front  and  back  mirror  layers.  The  sample  was  grown  at  the  University  of 
Arizona  in  June  of  1991.  Locations  of  lasing  and  non-lasing  spots  measured  are  contained 
in  figure  5.7  below. 


Figure  5.7.  Sample  2:  3  A  Active  Region  VCSEL  Geometry 


5.3.1  Material  Parameters  The  material  parameters  for  this  sample  are  sum¬ 
marized  below: 


Number  of  Layers  (front): 

44 

Number  of  Layers  (back): 

51 

Active  Region: 

GaAs  UGS 

Active  region  thickness 

.726  HTd 

Substrate: 

GaAs 

Design  wavelength: 

875  nm 

Low  Index  Layers: 

AlAs 

Calculated  Index  ®  875  nm: 

2.99573 

Physical  Layer  width: 

73.2  nm 

High  Index  Layers: 

Al.i37Ga.g73As 

Calculated  Index  <9  875  nm: 

3.52773 

Physical  Layer  width: 

61.8  nm 

The  calculated  spectral  reflectivity  of  Sample  2  based  on  the  parameters  described 
above  (with  no  loss  assumed)  is  represented  in  figure  5.8.  The  minimum  nearest  the  stop 
band  is  at  fa  824nm 
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Figure  5.8.  Sample  2:  Theoretical  Spectral  Reflectivity  near  875  nm 

5.S.2  Sample  2  Measurements  Figure  5.7  shows  the  spots  on  sample  2  which 
evidenced  lasing  and  non-lasing  behavior.  Again,  reflectivity  measurements  were  accom¬ 
plished  for  all  of  these  spots,  but  only  those  necessary  for  the  analysis  are  included. 

Sample  2;  Spot  3  : 

Lasing  Wavelength  885.8  nm 
Pump  Wavelength  840.4  nm 
Polarization:  Vertical 
Mode:  TEMio 

Figure  5.9  is  a  reproduced  image  of  the  TEMio  mode  this  spot  produced.  The  spot 
width  is  approximately  2.5  pm. 


Single  2  Thooreticftl  Reflectivity 
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VCSEL  Spot 


Figure  5.9.  Sample  2,  Spot  3:  TEMio  Lasing  and  Pump  Spots 


Reflectivity  As  for  sample  1,  the  reflectivities  measured  for  sample 
2,  spot  3  indicate  that  the  lasing  and  neighboring  non-lasing  areas  on  the  wafer  surface 
have  very  nearly  the  same  reflectivities.  Reflectivities  for  Sample  2  were  collected  every 
2.5  nm  to  ensure  minima  in  the  curve  were  not  bypassed. 


S«npl«  2Spai  3:  Refleclivitiet  for  Lasing  and  Non-Lasing  Spott 


Figure  5.10.  Sample  2  Measured  Spectral  Reflectivities  for  Lasing  and  Non-Lasing  Loca¬ 
tions  near  Spot  3 
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Sample  2;  Spot  4 


Lasing  Wavelength  879.0  nm 
Pump  Wavelength  833.3  nm 
Polarization:  Vertical 
Mode:  TEMoo 


Figure  5.11  shows  the  TEMoo  mode  this  spot  produced.  The  spot  width  is  approxi¬ 
mately  1.5  pm.  VCSEL  cipot 


Figure  5.11.  Sample  2,  Spot  4:  TEMoo  Lasing  and  Pump  Spots 


Reflectivity  The  spectral  reflectivity  for  Spot  4  shown  in  figure 
5.12  verifies  the  close  reflectivity  match  shown  for  spot  3  of  sample  2  above. 


Sample  2;  Spots  1  and  2  : 


Spotl: 


Lasing  Wavelength  882.7  nm 
Pump  Wavelength  838.2  nm 
Polarization:  Vertical 
Mode:  TEMoo 
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Figure  5.12.  Sample  2  Measured  Spectral  Reflectivities  for  Lasing  and  Non-Lcising  Loca¬ 
tions  near  Spot  4 


Spot  2: 

Lasing  Wavelength  880.5  nm 
Pump  Wavelength  835.8  nm 
Polarization;  Vertical 
Mode:  TEMqo 

Sample  2;  Non-Lasing  Spots  :  Spots  5  and  6  in  figure  5.7  did  not  lase. 
Reflectivities  for  these  spots  are  included  as  representative  of  the  non-lasing  areas  of  '•ample 
2.  The  reflectivity  data  for  these  areas  is  contained  in  figure  5.13. 

Sample2:  Reflectivity  Measurements  at  Non-Oblique  Incidence  The  mea¬ 
surements  presented  so  far  were  performed  with  the  focusing  lens  in  place  to  allow  reflec¬ 
tivity  measurements  of  small  spots  on  the  surface  of  the  wafers.  To  allow  comparison  to 
the  theoretical  reflectivity  for  normal  incidence,  the  focusing  lens  was  removed  and  the 
pump  beam  allowed  to  strike  the  si’rface  of  the  wafer  without  interference.  The  resultant 
larger  (~  3mm)  reflected  spot  forms  an  average  reflectance  over  a  larger  area  of  the  wafer. 
Two  measurements  were  performed  on  Sample  2  in  this  configuration;  one  at  an  area  on 
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Figure  5.13.  Sample  2  Measured  Spectral  Reflectivity  for  Non-Lasing  Locations:  Spots  5 
and  6 

the  wafer  with  many  lasers,  and  another  at  an  area  with  no  lasers.  This  measurement  was 
not  possible  for  sample  1,  which  was  much  smaller  than  sample  2,  and  scattered  much  of 
the  3  mm  pump  spot.  The  reflectivities  are  contained  in  figure  5.14 

5.3.3  Sample  2  Analysis  The  pump  and  lasing  wavelengths  for  all  four  spots 
indicate  higher  wavelengths  than  predicted.  The  indications  are  that  the  actual  lasing 
wavelength  for  the  structure  is  nearer  885  nm  than  the  designed  875  nm. 

The  reflectivity  curves  for  spots  3  and  4  indicate  a  close  match  in  reflectivities  between 
adjacent  lasing  and  non-lasing  spots.  The  quality  of  the  sample  2  wafer  was  appreciably 
better  than  sample  1,  which  was  an  earlier  fabrication,  has  fewer  mirror  layers,  and  has 
been  subjected  to  more  handling.  These  and  the  non-oblique  incidence  spectral  reflectivity 
curve  indicate  higher  pump  wavelengths  than  predicted  for  the  structure.  A  plot  of  the 
non-oblique  incidence  measurements  and  a  theoretical  spectral  reflectivity  curve  for  an  887 
nm  GaAs  VCSEL  is  shown  in  figure  5.15.  The  closer  match  is  evident. 

The  non-oblique  incidence  (measurements  with  no  focusing  lens)  spectral  reflectivi¬ 
ties  indicate  the  non-lasing  area  has  a  significantly  reduced  reflectivity  about  the  desired 
lasing  wavelength.  However,  the  reflectivities  obtained  with  the  focusing  lens  show  that 
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Figure  5.14. 


Figure  5.15. 
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Sample  2  Measured  vs  Theoretical  (875  nm  VCSEL)  Spectral  Reflectivities 
for  Lasing  and  Non-Lasing  Areas 
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Sample  2  Measured  Spectral  Reflectivities  for  Lasing  and  Non-Lasing  Areas 
vs  Theoretical  (887  nm  GaAs  VCSEL) 
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reflectivities  foi  adjacent  lasing  and  nonlasing  spots  are  nearly  the  same.  Spectral  reflec¬ 
tivity  measurements,  then,  can  give  a  general  idea  of  the  lasing  prospects  for  the  wafer, 
but  do  not  guarantee  a  specific  spot  will  lase. 

All  polarizations  were  linear  and  indicated  a  90®  rotation  to  those  measured  for 
sample  1.  Again,  no  information  regarding  the  growth  direction  was  available,  and  it  is 
likely  this  sample  provides  polarization  orthogonal  to  that  in  sample  1  because  the  same 
physical  mechanisms  apply,  but  the  sample  is  rotated  90®  on  the  measurement  stand. 

As  noted  in  the  spot  descriptions,  both  TEMio  and  TEMqo  modes  were  found  on  the 
sample.  Pictures  of  these  modes  are  provided  in  figure  5.9  and  5.11  above.  The  TEMio 
lasing  mode  occurred  on  spot  3.  The  mode  was  stable,  remaining  unaltered  for  changes  in 
pump  intensity,  wavelength,  and  position. 
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5.4  Design  VCSEL  Back  Mirror 


The  initial  back  mirror  was  designed  to  measure  the  quality  of  the  MBE  fabrica¬ 
tion  before  proceeding  to  the  full  VCSEL  structure.  The  desired  lasing  wavelength  and 
reflectivity  between  the  GaAs  spacer  and  the  back  mirror  were  the  primary  concern.  The 
plot  for  the  reflectivity  in  the  GaAs  spacer  is  in  figure  3.4  on  page  3-8.  For  experimental 
measurements,  the  back  mirror  reflectivity  was  to  be  determined  with  the  incident  medium 
being  air,  rather  than  GaAs.  To  stave  off  oxidation  of  the  AlAs  layers,  a  top  A/4  layer 
was  deposited  on  the  mirror  stack,  resulting  a  a  17  period  stack  of  AlAs/GaAs  pairs  on  a 
GaAs  substrate. 


5.4-1  Material  Parameters  The  parameters  for  the  bottom  mirror  in  air  were: 

Number  of  Layers:  34 

Substrate:  GaAs 

Peak  Reflectivity:  99.147  %  @  920  nm  in  air 


Low  Index  Layers:  AlAs 
Calculated  Index  ®  920  nm:  2.98151 
Physical  Layer  width:  77.1  nm 


High  Index  Layers: 
Calculated  Index  ®  920  nm: 
Physical  Layer  width: 


GaAs 
3.5727 
64.3  nm 


Wafer  diameter:  50  mm 


5.4.2  Back  Mirror  Measurements  The  spectral  reflectivity  curve  in  air  for  such 
a  structure  was  computed  and  the  curve  is  plotted  in  figure  5.17  along  with  the  measured 
data  for  the  grown  wafer.  The  reflectivity  attained  by  the  wafer  reaches  the  desired  re¬ 
flectivity  for  the  bottom  mirror.  Since  the  calculated  spectral  reflectivity  curves  do  not 
include  absorption  and  scattering  losses,  the  reflectivity  obtained  by  the  wafer,  can’t  be 
expected  to  reach  the  predicted  curve.  The  degree  to  which  the  shape  of  the  measured 
spectral  reflectivity  curve  matches  the  design  curve  indicates  the  proximity  of  the  layer 
thicknesses  to  design  thicknesses.  The  stop  band  width  of  the  mecisured  curve  can  indicate 
if  the  indices  of  refraction  of  the  layers  are  close  to  those  specified,  as  discussed  in  chapter 
II. 
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Initially,  naeasurements  were  taken  at  two  locations  on  the  wafer;  one  set  of  data  at 
the  wafer  center,  and  another  at  10  mm  radially  out  from  the  center.  The  configuration 
in  figure  4.3  was  used  for  these  measurements,  as  that  configuration,  with  the  broad  Ti- 
Sapphire  wavelength  range,  provided  a  more  complete  spectral  reflectivity  measurement. 

Wafer  Center  Spectral  Reflectivity  :  The  measurements  shown  in  figure 
5.16  were  taken  at  a  7®  oblique  incidence,  which  shifts  the  spectral  reflectivity  curve  to 
slightly  shorter  wavelengths  than  for  normal  incidence  («  5  nm).  The  first  minimum  next 
to  the  stop  band  is  very  close  to  the  desired  wavelength,  and  the  stop  band  width  is  nearly 
the  same  as  predicted.  This  indicates  the  layer  thicknesses  and  indices  of  refraction  are  as 
needed  at  the  center  of  the  wafer. 


700  730  800  SSO  900  9S0  lOOO 

wavcleniUi  (tan) 


Figure  5.16.  Spectral  Reflectivity  Measurements  at  Center  of  Designed  920  nm  Back 
Mirror 


Wafer  10  mm  Spectral  Reflectivity  :  The  measurements  at  a  location  10 
mm  out  from  the  wafer  center  were  taken  at  a  65®  oblique  incidence.  This  high  angle  of 
incidence  shifts  the  spectral  reflectivity  curve  (w  28  nm  )  down.  The  spectral  reflectivity 
curves  here  indicate  a  close  match  for  the  edge  of  the  stop  band  and  reproduce  the  secondary 
maximum  at  970  nm  well.  No  conclusions  about  the  indices  of  refraction  can  be  drawn. 
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since  this  curve  has  not  captured  both  edges  of  the  stop  band,  but  it  is  unlikely  the  index 
for  each  layer  will  change  with  location  on  the  wafer  and  the  stop  band  in  figure  5.16 
indicates  the  designed  refractive  indices  are  present. 


Detigned  Back  Minor  Meaaured  vt  Theoretical  Rdlectivity  at  65  degreea  OUique  Incidence 


Figure  5.17,  Spectral  Reflectivity  Measurements  at  10  mm  radial  out  from  Center  of 
Designed  920  nm  Back  Mirror 


Wafer  Radial  Measurement  Comparisons  :  Comparison  of  reflectivities 
at  locations  on  the  wafer  radially  out  from  the  center  were  next  completed.  Figure  5.18 
contains  a  comparison  of  reflectivities  for  the  center  of  the  wafer,  and  for  locations  at  10 
mm,  and  20  mm  radially  out.  Figure  5.19  is  a  comparison  of  two  spots  located  5  mm  out 
from  the  center  of  the  wafer.  This  is  an  indication  of  the  uniformity  of  the  growth  at  a 
constant  radius. 

Figure  5.18  shows  a  large  variance  in  the  reflectivity  of  the  wafer  as  the  spot  is  moved 
from  the  center  of  the  wafer  radially  out  to  points  10  and  20  mm  from  center.  The  stop 
band  tends  to  move  to  shorter  wavelengths  indicating  the  layers  are  becoming  thinner. 
The  high  edge  of  the  stop  band  near  the  edge  (20  mm  from  center  of  wafer)  has  moved  to 
940  nm.  This  is  approximately  35  nm  from  the  calculated  location  of  40%  reflectivity  for 
a  7®  oblique  incidence  (see  Calculated  curve  in  figure  5.16).  If  this  shift  is  typical  for  this 
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radial,  the  Bragg  frequency  has  moved  to  885  nm,  indicating  the  mirror  layer  thicknesses 
are  reduced. 


S40  S60  Sto  900  920  940  960 

«nv<lcn(th  (am) 


Figure  5.18.  Spectral  Reflectivity  Measurements  at  Center,  10mm  radial,  and  20  mm 
radial  of  Designed  920  nm  Back  Mirror 


To  determine  if  the  mirror  layers  were  uniform  on  a  constant  radius  from  the  center, 
two  points  on  the  wafer  5  mm  from  the  center  in  opposing  directions  were  measured. 
Uniform  thickness  on  a  constant  radius  is  expected,  since  the  wafer  is  rotated  during  the 
MBE  growth  specifically  to  improve  this  constant  radius  uniformity.  The  result  of  this 
measurement  is  shown  in  figure  5.19,  and  indicates  the  reflectivities  for  the  two  locations 
are  very  nearly  the  same.  This  radial  uniformity  can  help  to  enlarge  the  areas  of  the  wafer 
having  the  desired  layer  thicknesses 
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Refleaivity 


Coinpwuaa  of  Ataaolute  Refloctivitiei  for  Back  Mkror  of  AlAa/CaAi  (34  Uyen):  Sample  #34S  (two  radiala) 


MO  KO  UO  900  920  940  960 

wavelanfth  (an) 


Figure  5.19.  Spectral  Reflectivity  Measurements  at  two  radials,  5mm  from  Center  of 
Designed  920  nm  Back  Mirror 
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VI.  Conclusion 


The  pace  of  improvements  to  the  basic  VCSEL  design  is  continuously  accelerating. 
This  research  verified  the  capability  to  locally  produce  VCSEL-quality  structures  and  char¬ 
acterize  their  performance.  The  by-product  of  this  effort  is  an  initial  in-house  capability  for 
design  and  characterization  of  VCSEL  structures;  and  a  nearby  means  for  fabrication.  The 
research  purposely  centered  on  optically  pumped  devices,  due  to  the  means  of  production 
available  locally  and  the  structures  we  could  “borrow”  from  the  University  of  Arizona. 

6.1  Characterization 

The  laboratory  configurations  used  for  device  characterization  are  detailed  in  chapter 
IV.  The  configuration  in  figure  4.1  provided  the  means  to  optically  pump  the  VCSELs 
supplied  by  the  University  of  Arizona,  measure  output  beam  characteristics,  and  measure 
spectral  reflectivity  from  the  wafer.  The  relatively  narrow  spectral  range  of  the  dye  laser 
prompted  additional  measurements  with  the  Ti-Sapphire  system  in  figure  4.3.  The  Ti- 
Sapphire  system  provided  a  means  to  cover  the  entire  stop  band  of  the  Fabry-Perot  curve 
for  the  VCSELs. 

Measurements  were  accomplished  with  and  without  a  focusing  lens.  The  lens  was 
necessary  to  concentrate  the  pump  on  the  VCSELs.  The  focusing  lens  also  provided  a 
means  to  isolate  the  reflectivity  from  an  individual  VCSEL  structure.  The  convolving  effect 
of  the  lens  on  the  spectral  reflectivity,  however,  was  unwanted.  Measurements  without  the 
focusing  lens  provided  the  spectral  reflectivity  for  a  broad  area  (3mm)  relative  to  the 
VCSEL  diameter  (~  3  -  20  /x  m)  without  the  convolving  effect  on  the  spectral  reflectivity. 
These  measurements,  however,  average  the  reflectivity  of  lasing  and  non-lasing  spots  on 
the  wafer.  Conclusions  drawn  and  areas  that  warrant  further  investigation  follow. 

6.1.1  Spectral  Reflectivity  Reflectivity  measurements  from  lasing  and  non-lasing  ar¬ 
eas  of  the  samples  indicate  a  cursory  spectral  reflectivity  measurement  can  quickly  exclude 
samples  with  very  poor  reflectivity  (<  80  %).  The  measurements  with  focused  pump  spot 
show  both  lasing  and  non-lasing  locations  on  the  wafer  can  have  nearly  the  same  spectral 
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reflectivity.  The  reason  why  one  area  lases  while  an  adjacent  area  does  not  should  be 
investigated  further.  There  may  be  multiple  mechanisms  at  work. 

Surface  Roughness  The  scattering  from  the  surface  of  the  wafer  may 
contribute  to  this  effect,  scattering  more  of  the  incident  pump  at  the  non-lasing  location. 
One  would  expect  a  lower  reflectivity  from  the  area  with  more  surface  scatter,  since  not  all 
the  scattered  light  can  be  expected  to  reflect  back  toward  the  detector.  Here,  the  action  of 
the  focusing  lens  (see  figure  4.2)  may  act  to  collect  a  large  part  of  the  scattered  light  and 
provide  no  indication  of  how  much  of  the  light  incident  at  the  surface  is  actually  entering 
the  active  region.  When  the  pump  intensity  is  increased,  the  scattering  still  results  in  a 
reduced  pump  intensity  in  the  active  region,  possibly  insufficient  for  the  onset  of  lasing. 
A  closer  inspection  of  the  surface  roughness  of  the  sample,  with  regard  to  lasing  and  non¬ 
lasing  areas,  would  be  in  order.  Future  studies  will  need  a  precise,  repeatable  means  of 
determining  position  on  the  wafer  as  the  sample  is  moved  from  the  lasing  configuration 
to  the  microscopic  analysis  configuration  to  explore  this  phenomenon.  The  laboratory 
configuration  in  figure  4.1  provides  precision  and  repeatability,  but  only  if  the  sample  is 
not  moved  from  the  positioner.  This  configuration  does  not  lend  itself  to  microscopic 
viewing  of  the  sample. 

Guide  Mechanisms  Another  mechanism  to  explore  is  the  gain/index 
guiding  used  to  stimulate  these  optically  pumped  VCSELs.  The  use  of  injected  carri¬ 
ers  to  self  guide  the  lasing  process  may  be  affected  by  surface  roughness  also.  Here,  the 
pump  input  may  be  generating  a  sufficient  amount  of  carriers  in  the  active  region,  but  the 
proximity  of  the  active  region  to  a  more  favorable  exit  spot  may  tend  to  channel  the  beam 
to  that  spot  (see  figure  6.1) 

This  prompts  a  reason  for  a  phenomenon  noted  in  the  lab  experiment.  The  imaged 
output  of  the  laser  spot  tended  to  shift  slightly  {ss  l^m)  as  the  pump  beam  passed  across 
the  lasing  location  on  the  wafer.  As  the  spot  moved  off  the  optical  axis,  the  intensity 
decreased  until  lasing  ceased,  leaving  a  photoluminescent  (but  non-lasing)  spot  shifted 
back  to  the  optical  axis. 
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lavcmd  Po^poUtioB 


Figure  6.1.  Non-Coincident  Pump  and  Exit  Beam  Locations  due  to  Surface  Roughness 


If  surface  roughness  provides  a  more  favorable  exit  point  adjacent  to  the  pump  spot, 
the  exiting  beam  would  move  to  that  location,  feeding  off  the  carriers  generated  at  the 
pump  location.  With  the  lab  configuration  used  for  this  effort,  the  focusing  lens  would 
image  the  exit  beam  at  a  point  removed  from  the  optical  axis,  as  shown.  The  further  the 
pump  moves  from  this  more  favorable  exit  point,  the  more  loss  encountered  by  the  carriers 
(they  are  on  the  fringes  of  the  high  gain  region,  and  experience  more  loss),  eventually, 
encountering  too  much  loss  to  support  lasing. 

The  ability  to  image  larger  areas  on  the  wafer  would  enhance  this  study.  The  ir¬ 
regularities  (especially  the  “islands”  mentioned  in  chapter  5  in  the  spot  2  analysis)  were 
approximately  the  same  size  as  the  pump  spot,  making  it  difficult  to  compare  the  irregular¬ 
ities  to  normal  structures  nearby.  These  “islands”  were  also  witnessed  in  VCSEL  imaging 
experiments  at  the  University  of  Arizona  [17],  and,  since  they  tend  to  exhibit  lasing,  may 
provide  useful  information  for  future  growths. 

6.1.2  Polarization  With  complete  wafer  growths  available  in  the  future,  another 
possible  avenue  of  investigation  is  the  link  of  polarization  and  growth  patterns.  As  dis¬ 
cussed  in  chapter  II  the  link  between  strained  active  regions  and  polarization  has  been 
established.  The  close  lattice  match  of  GaAs  to  the  Al*Gai_,As  ternaries  provides  very 
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little  strain.  In  his  investigations,  Jewell  [21]  found  apparently  random  orientation  to  the 
linear  polarization  output  beams  from  an  optically  pumped  GaAs  VCSEL.  Further  in¬ 
vestigation  into  this  area  might  show  a  link  between  polarization  direction  and  growth 
patterns. 

6.1.3  Laboratory  Configuration  Improvements  Many  of  the  improvements 
already  planned  will  significantly  improve  the  capability  for  characterization  and  VCSEL 
stimulation. 

-  -  The  dye  laser  used  in  this  experiment  was  difficult  to  maintain,  had  a  limited 

frequency  range,  and  exhibited  output  power  fluctuations  which  made  it  dif¬ 
ficult  to  capture  reliable  reflectivity  data.  The  use  of  the  Ti-Sapphire  system 
highlighted  these  shortcomings.  The  seven  measurements  used  in  this  report 
from  the  configuration  with  the  Ti-Sapphire  laser  (see  figure  4.3)  were  done  in 
one  day.  The  dye  laser  configuration  would  have  taken  at  least  five  days  to 
accomplish  the  same  number  of  me<tsurements  and  would  not  have  been  able  to 
cover  the  desired  frequency  range. 

-  Faist  [14]  suggests  transmission  electron  microscopy  (TEM)  is  essential  for  more 
detailed  analysis  of  the  surfaces  involved.  Any  future  study  into  the  surface 
roughness  aspects  of  the  problem  will  rely  greatly  on  the  availability  of  some 
form  of  microscopy  for  physical  comparison  to  the  lasing  structure’s  behavior. 

-  The  addition  of  an  optical  multichannel  analyzer  will  significantly  speed  the 
spectral  reflectivity  measurements  by  using  a  broadband  light  source  rather 
than  a  laser  to  perform  these  measurements. 

6.2  Design 

Design  improvements  can  be  accomplished  in  many  areas.  The  vertical  cavity  struc¬ 
tures  described  herein  are  intentionally  entry  level.  Some  immediate  improvements  in  the 
efficiency  of  the  VCSELs  grown  can  be  realized,  some  are  longer  term. 


6-4 


1.  Incorporate  separate  confinement  heterostructure  (SCH)  into  the  quantum  well  de¬ 
sign  described.  This  will  improve  the  longitudinal  confinement  factor  (F)  for  the 
active  region,  and  would  involve  only  a  few  extra  growth  layers. 

2.  Graded  index  (GRIN)  quantum  well  structures  to  improve  the  capture  of  the  carriers 
in  the  InGaAs  wells. 

3.  Higher  reflectivities  in  the  mirrors  should  be  attempted.  Higher  reflectivities  will 
favor  lower  thresholds  as  long  as  the  losses  in  the  mirrors  can  be  kept  low.  Coldren 
indicates  the  point  of  diminishing  returns  for  three  and  two  quantum  well  structures 
is  near  99.7  %  mean  mirror  reflectivity  [9]  and  goes  to  lower  reflectivities  for  more 
quantum  wells. 

4.  The  literature  is  very  sparse  on  the  material  properties  of  In^Gai.^As  on  GaAs 
substrates.  The  lack  of  this  data  inhibits  the  ability  of  the  designer  to  determine 
the  cavity  gain  and  loss  profiles.  Efficiency  of  the  design  depends  on  these  param¬ 
eters.  Many  of  the  major  successes  in  the  VCSEL  efficiency  and  power  areas  have 
resulted  from  an  in-house  capability  to  measure  prospective  material  parameters. 
Some  suggest  the  construction  of  simpler  in-plane  lasers  and  subsequent  parameter 
measurement  in  these  devices  to  provide  the  necessary  parameters  for  an  efficient 
VCSEL  design.  This  approach  may  be  an  appropriate  meaisure  if  efficiency  is  to  be 
pursued. 

5.  Ion  implantation  for  gain  guiding  and  etching  for  index  guiding  should  be  attempted. 
A  comparison  to  an  identically  grown  VCSEL  wafer  without  etching  or  ion  implanta¬ 
tion  could  indicate  if  these  processes  provide  an  essential  extra  carrier  confinement, 
and  increase  the  number  of  usable  lasers  on  the  wafer. 

6.  The  incorporation  of  electrical  stimulation  of  the  active  medium,  two  dimensional 
arrays,  coupled  outputs,  and  visible  wavelength  structures  make  VCSEL  fabrication 
an  area  with  an  extremely  big  potential  return. 
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Appendix  A.  Spectral  Reflectivity  FORTRAN  Programs 


The  programs  below  were  first  developed  at  the  University  of  Arizona,  and  were 
modified  to  apply  to  the  VCSEL  structures  and  DBR  mirrors  measured  for  this  thesis. 


A.l  Program  to  Calculate  the  Spectral  Reflectivity  in  Fabry-Perot  Etalon 
with  Distributed  Bragg  Reflector  Mirrors  (VCSEL  Structure) 
c  program  vcsel.l 

c  This  program  calculates  the  transmission  and  reflection  of  a  GaAlAs 
c  Fabry-Perot  etalon  sith  integrated  mirrors  as  a  function  of  the  input 
c  beam  vavelength. 

PROGRAM  vcsel 

real  na,ns,imbdal,lmbda2.1mb 

real  lmbda,TRARS(1000) ,cosd(100) .Imbdal ,lmbda2 

real  Imbdac,  indxcL,  indxcH,  indxcs ,indxH,  indxL,indxs 

character's ,  thick , oblique , pol , reflet , yorn , meix 

complex  indx(lOO) ,phase(100) ,cindxh,cindxl,cindxs 

complex  r ,t ,s ,u,rr ,tt ,ss ,uu,rrr ,ttt,ss3,uuu,b,c 

common  sina.lmbda 

c  This  section  initializes  the  parameters  for  the  menu, 

c  allowing  the  star-u^  program  to  have  a  set  of  parameters 

c  which  will  produce  a  default  set  of  data, 

c 

WRITE(*,2) 

2  FORMATC’  ',’The  reflectivity  data  is  in  FP.dat’) 

write(e,4')  ’The  current  default  parameters  2u:e: ' 
c  Humber  of  front  and  back  mirror  layers 

ILF=22 
RLR=33 

c  Spacer  length 

spcerl=.515 

c  Loss  in  the  layers 

al=0.0 

c  X  concentrations  in  the  low,  high  index  layers  emd  the  spacer 

xl=1.0 
xh=0.0 
xs=0.0 

c  Whether  to  allow  changes  in  the  layer  thicknesses, 

c  oblique  incidence(which  requires  a  pol^u^ization) 

thick=’no’ 
oblique=’no’ 
pol=’TE’ 

c  Lasing  wavelength 


A-: 


lmbdac= . 92 

c  Beginning,  ending,  and  spacing  wavelength  lor  plot 

lmbdal=.7 
lmbda2=l. 
nlinbda=300 

c  Default  layer  Thickness  variation  factors,  high  index  layers, 

c  losindex  layers,  and  spacer 

lctorl=l . 
fctorh=l . 
f ctors=l . 

c  The  incident  medium  for  the  reflection 

reflct='air’ 

c  Whether  to  remain  in  the  program 

yom=’y’ 

c  This  part  of  the  program  produces  the  menu  on  the  screen 

10  open(3,  file='FP.dat') 

write(6,20)  ’[1]  Ho.  of  layers  front', ’[2]  Ho.  of  layers  back' 
Brite(6,25)  HLF,HLR 
20  F0RMAT(10X,A,  lOX,  A) 

25  FORMAT (14X, I,  16X,  I) 

print*,  '  ' 

write(6,30)  'C3]A1  concentrations; ', 'high  n  layer', 

*  'low  n  layer',  'spacer  layer' 

vrite(6,35)  xh,zl,xs 
30  F0RNAT(1X,A,  6X,  A,5X,A,5x,A) 

35  F0RMAT(30X,f6.2,  12X,  15.2, llx,f5.2) 

print*,  '  ' 

write(6,40)  '[4]  Design  Havelength(nm) ' , 
k  ' [5]  Spacer  Thickness  ' 

write(6,42)  lmbdac*1000,spcerl 
40  F0RMAT(2X,A,  lOX,  A) 

42  FaRMAT(12X,f6.2,  30X,  16.3) 

writeC*,*) 

write(6,4S)  '[6]  Oblique  Incidence', 
k  'Polarization  ' 

write (6, 47)  oblique,  pol 
45  F0RNAT(2X,A,  15X,  A) 

47  F0RMAT(12X,A,  27X,  A) 

write(*,*) 
write(*,*) 

writa(6,50)  '[7]  Plot  Data  in  FP.dat' , 'start  (nm) ' , 
k  'end  (nm)',  'number  of  points' 

write(6 ,55)  lmbdal*1000 ,lmbda2*1000 ,nlmbda 
50  F0RMAT(1X,A,  5X,  A,2X,A,2x,A) 

65  FORMAT (29X, 17.1,  4X,  f7.1,2x,I) 

write(*,*) 
write(*,*) 
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vrita(6,60)  ’[8]  Vary  layer  thickness' 

Brite(6,65)  thick 
60  FORMAT (IX. A) 

66  FORNATClOX.A) 

Brite(e,*) 

Brite(*,*) 

Brite(6,67)  ’[9]  Absorption  constant' 

Brite(6,68)  al 

67  FORMATdX.A) 

68  FORMAT (lOX, 16. 2) 

BriteC*,*) 

Brite(*,*) 

c 

c  Here  the  interaction  Bith  the  user  starts,  allouing  changes 

c  to  parameters  on  the  menu 

c 

printe,  'Type  in  the  parameter  to  change  (one  parameter 
k  at  a  time,  0  to  proceed' 
read(*,>»)choice 
il ( choice . eq . 1 ) then 

Brite(*,*)  'number  ol  layers  ol  the  mirror=?‘ 
read(*,*)  nil 
endil 

il(choice . eq. 2)then 

Brite(*,*)  'number  ol  layers  ol  the  mirror=?' 
read(*,*)  nlr 
endil 

il (choice . eq. 3)then 

Brite(*,s)  'Al  content  ol  high-index  layer=?' 
read(*,*)  xh 

BriteC*,*)  ’Al  content  ol  loB-index  layer=?’ 
read(*,*)  xl 

Brite(*,*)  'Al  content  ol  spacer' 
read(*,*)  xa 
endil 

il (choice . eq.4)then 

Brite(*,*)  'design  Bavelength  ol  mirror  in  nm=?’ 
read(4',e)  Imbdac 
lmb<lac=lmbdac* .  001 
endil 

il (choice. eq. 5) then 

Brite(*,*)  'spacer  thickness=?' 
read(*,*)  spcerl 
endil 

il ( choice . eq . 6 . and . oblique . eq . ' no ' ) then 
oblique='ye8' 

Brite(*,*)  'Incidence  angle  in  degrees?' 
read(*,*)  angle 

Brite(*,*)  'Polarization,  il  TE,  input  1:' 
read(*,*)  npol 
il  (npol .  eq .  0}  then 
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pol='TM' 

else 

pol='TE' 

endil 

elseiK choice. eq. 6. and. oblique. eq. 'yes'} then 
obliques 'no' 
anglesO . 
npol=l 
pols'TE' 
endif 

il (choice . eq. 7) then 

vriteC*,*)  'Imbdal  (beginning  wavelength:  nm)' 

read(*.*)  Imbdal 

Imbdalslfflbdale.OOl 

write(*,*}  'Imbdal  (ending  wavelength:  nm) ' 
read(*,v)  lmbda2 
Imbda2=lmbda2v . 001 

write(v,v)  'nlmbda  (number  of  points)' 
read(e,e)  nlmbda 
endif 

if (choice . eq. 8 . and . thick. eq. 'no ' )then 
thicks 'yes' 

write(v,*)  'Thickness  multiplier  for  High  index  layers?' 
read(e,*)  fctorh 

write(*,*)  'Thickness  multiplier  for  Low  index  layer*?' 
read(*,*)  fctorl 

write(*.*)  'Thickness  multiplier  for  spacer*?' 
read(v,*)  fetors 

elseif (choice . eq . 8 . and . thick . eq . 'yes ' ) then 
thick* 'no' 
fctorl* 1 . 
fctorh*! . 
fetors*!, 
endif 

if ( choice . eq . 9 . and . reflet . eq. ' cavity ' )then 
write(*,*)  'absorption  constant*?' 
road(*,*)  al 
endif 

if (choice. ne.0)goto  10 
c 

c  initializing  constants  for  this  run 

c 

PI*3. 141592 
angle*angleepi/ 180 . 
sina*sin(angle) 
na=l . 

conv*6 . 626e-34*3 . el4/l . 6e-19 

c  determine  the  width  of  high- index  (indxH)  layer  (widthH)  and 
c  low-index  (indxL)  layer  (widthL) .  lamdac  is  the  design  frequency  of 
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c  the  etalon. 


hvc^conv/lmbdac 

CALL  REF(xl.  hvc,  indxcL) 

CALL  REF(xh,  hvc.  indxcH) 

CALL  REF(XS,flVC.IHDXCS) 

SPCERL=SPCEIlL*FCTORS 
vidthH  s  lmbdac/4 . /indxcHepCTORB 
vidthL  =  lmbdac/4./indxcLeFCT0RL 
print*, 'spacer  index  =  indxcs 
WRITE(*,*)  'change  mirror  thickness?  il  yes,  input  1' 

READ(*,*)  number 
i2 (number. eq.l)  then 

WRITE(*.*)  'THICKHESS  oi  HIGH  INDEX  LAYER=?' 

READ(*,'(F5.3)')  widthh 

HRITE(*.*)  'THICKHESS  oi  LOW  INDEX  LAYER=?' 

READ(*. '(F5.3)')  widthl 
endif 

c  prints  sidth  of  high,  los  and  spacer  layers  for  design 

«rite(*,lll)  vidthh.uidthl.SPCERL 

111  format (lx. 'high  in  um=’ ,f7.6,5x, 'los=* ,f7.6,Sx. 'spacer=' ,f8.6) 
WRITE(4,76)  tfIDTHH, WIDTHL 
WRITE(4,80)  INDXCH, INDXCL. INDXCS 
75  FORMAT(//, 

+1X,F5.3,BX,'HIGHINDX  THCKNSS  ' ,F6.4,/, 'LOW  INDX  THCKNSS  ',F6.4) 

80  FORMAT(//,'HIGH  INDX  INDEX  ' .F6.3,5X, 'LOW  INDX  INDX  '.F5.3,/, 

+' SPACER  INDEX  '.F5.3) 

c  calculate  the  transmission  trans  of  the  etalon  as  the  wavelength  is 
c  varied  from  Imbdal(um)  to  lmbda2(um) . 

XX1=0. 

XX2=0. 

do  200  ksl.nlmbda-t-l 

lmbda=lmbdal * (lmbda2-lmbdal ) /float (nlmbda) *f loat (K- 1 ) 
lmb=lmbda* 1000 

c  alphal  is  the  absorption  constant  multiplied  by  the  thicknes.s. 
alphal=0 . 

if(  Imbda  .It.  0.875  )  alphal=al 

plmbdaslmbda/4/pi 

hv=conv/lmbda 

CALL  REF(xl,  hv.  indxL) 

CALL  REF(xh.  hv,  indxH) 
call  ref(xs,  hv,  indxs) 

CALL  REF(0.,HV,NS) 
xxx=-alphal*lmbda/4/pi 
cindxsscmplx ( indxs , zxx) 
xxl^xxxeO. 

cindxh*cmplx ( indxh , xx 1 ) 
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cindzl^cmplx ( indxl . 0 . ) 

c  call  tha  phase  change  phase(m)  by  passing  thru  each  layer  at 
c  wavelength  Imbda. 

do  201  m  =  2. nil, 2 

indz(m)  =  cindxL 

call  pha(indz(m) .widthL.cosd(m) ,phase(m}) 

201  continue 

do  202  m  =  l,nll,2 

indz(m)  =  cindxH 

call  pha(indx(m) .widthH,cosd(si) .phase (m)) 

202  continue 

m=nlf+l 

indx(m)=cindxs 

call  pha(indx(in) ,spcerl,cosd(m) ,phase(m)) 

do  203  m  =  nll-i-2,nll->-nlr+l ,2 
indx(m)  =  cindxL 

call  pha(indx(m) ,widthL,cosd(m) ,phase(m)) 

203  continue 

do  204  B  =  nll+3,nll+nlr+l,2 
indx(m)  ~  cindxH 

call  pha ( indx (m) , widthH , cosdCm) , phase (m) ) 

204  continue 

c  calculate  the  transmission  trans  at  wavelength  Imbda. 
c  refer  to  Thin-Film  Optics  (2/e}  by  A.  Nacleod. 

if (numbera.eq. 1)  then 
do  205  B=l,nlr-«-nlf-*-l 
if (npol.eq. 1)  indx(m)=indx(B}ecosd(m) 
if (npol.ne. 1)  indx(B}=indx(B)/cosd(B) 

205  continue 
endif 

r  =  ccos(phase(l)) 

s  =  (0,1.)  *  csin(phase(l))  /  indx(l) 

t  =  (0,1.)  *  csin(phase(l))  *  indx(l) 

u  =  r 

do  206  B=2,nlr-i-nlf'*-l 

rr  =  ccos(phase(a)) 

ss  =  (0,1.)  ♦  csin(phase(B))  /  indx(B) 
tt  =  (0,1.)  *  csin(phaae(B))  *  indx(m) 
uu  =  rr 
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206  continue 

b=r+s*n8 

c=t+u*na 

c  lacr  is  to  include  the  ellect  by  reflection  at  the  substrate-air 
c  interface. 

FACR= 1 . - ( (lA- 1 ) / (Il+l ) ) **2 

trans  (K)  =4 .  enaens/cabs  (naeb-t-c)  eeSePACR 

REFL=cab8 ( (naeb-c) / (na*b+c) ) **2 

c  Imbda  is  multiplied  by  1000  to  change  the  unit  from  um  to  nm. 
c  trans (k)  and  refl  are  multiplied  by  100  to  get  '!>  of  transmission. 
HRITE(3,260)  LMBDA*1000.,  REFL*100.,  tran8(k)*100. 

c  the  following  is  the  routine  to  find  peak  positions  and  the 
c  corresponding  peak  values  of  transmission. 

IF(K.GT.l)  XX1=TRAIS(K)-TRAHS(K-1) 
IF(XX2.GT.0.AID.XX1.L'..0.)  THE! 

WRITE(*.222)  PEAK.tRAHS(K-l) 

EID  IF 
XX2=XX1 
peak=LMB0A 
200  continue 

222  formatClz,  'peak  position  in  um=' ,f8.6,3x, 'transmission^' ,f8.6) 
2S0  F0RMAT(1X,3F10.5) 

CL0SE(3} 

writeC*,*) 
write(*,*) 
print* , ' continue? ' 
r*ad(*,*)yom 
if  (yom.eq.  'y')then 
goto  10 
endif 
EID 
c 

c  This  subroutine  calculates  the  phase  change  in  passing  thru 
c  each  layer  for  given  values  of  the  refractive  index  (  xindz  ) ,  the 
c  thickness  (  width  },  the  incident  angle, 
c 

subroutine  pha(xindx .width, cosd .phase) 
complex  xindx, phase 
common  sina.xlmbda 
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8inan=8ina/real (xindx) 
aa inansas in ( 8 inan ) 
cosd=co8 (as Inan) 

phase  =  2.  *  3.141692  *  xindx  *  sidth  *  cosd  /  xlmbda 
return 
end 
c 

c  This  subroutine  calculates  refractive  index  of  AlGaAs  for 
c  given  values  of  A1  composition  (  xconc  )  and  uavelength  (  ev,  in 
c  units  of  eV) . 
c 

SUBROUTINE  REF(Xconc.  ev.  XN) 

aO  -  6.3  -t-  19.0  e  xconc 
bO  =  9.4  -  10.2  *  xconc 

EO  =  1.426  +  1.165  ♦  xconc  +  0.37  *  xconc  *  xconc 
EOD  =  1.766  +  1.116  *  xconc  +  0.37  ♦  xconc  e  xconc 
xki  =  ev  /  EO 
xkiso  =  ev  /  EOD 

fx  =  (2.  -  sqrt(l.+xki)  -  sqrt(abs(l.-xki)))  /  xki  /  xki 

fxso  =  (2.  -  sqrtCl .-t’Xkiso)  -  sqrt(abs(l-xkiso)))  /  xkiso  /xkiso 

val  s  EO  /EOD 

epis  =a0*(fx+0.6e  val**1.5  ♦  fxso)  +  bO 
xn  =  sqrtCepis) 

XX  =  3.66  +  (ev  -  1.6)/. 3  ♦  .13 

if  ((ev.ge.l.4.or.xx.lt.xn).and.xconc.eq.O.)  xn  =  xx 

return 

end 


A. 2  Distributed  Bragg  Reflector  Mirror  Spectral  Reflectivity  Calculation  Pro¬ 
gram 

c  This  program  calculates  the  reflection  by  a 
c  GaAlAs  mirror  as  a  function  of  the  input  beam  uavelength. 

PROGRAM  DBR 
integer  nl, choice 

real  nspace,  nsub,  lambdal,  lafflbda2.top 
real  lambda,  TRAIS( 1000) ,cosd( 100) 
real  lambdac,  indexcL,  indexcH,  indexH,  indexL 
complex  index(lOO) ,phase(100) .cindexh.cindexl 
complex  r,t,8,u,rr,tt,ss,uu,rrr,ttt,sss,uuu,b,c 
charactered ,  borf , thick . oblique ,pol .reflet .yom .max 
common  sina, lambda 

c  Setting  menu  defaults 

c  Default  choice  of  back  mirror  to  be  the  calculated  one 

borf* 'back' 
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lumber  of  layers 

IL=33 

X  concentrations  in  tbe  lou,  high  index  layers  and  the  spscer 
xl=1.0 
xh=0.0 
xs=0.0 

Whether  to  allov  changes  in  the  layer  thicknesses, 
oblique  incidence (uhich  requires  a  polarization) 
thick='no* 
oblique=’no' 
pol='TE' 

Lasing  uavelength 
lambdac= . 92 

Beginning,  ending,  and  spacing  wavelength  for  plot 
lambdal=.7 
lafflbda2=l. 
nlambda=300 

Default  layer  Thickness  variation  factors,  high  index  layers, 
lowindex  layers,  and  spacer 
factorial. 
factorh=l . 

Incident  medium  for  the  reflected  light  (cavity  or  air) 
use  cavity  if  you  want  the  reflection  the  laser  sill  experience 
air  if  you  want  waht  to  expect  when  meas'ured  in  the  lab 
reflet* ’cavity’ 

Maximize  reflectivity,  will  automatiacally  throw  out  layers  that 
don’t  add  reflectivity  and  notify  you  it  did  so 
max* ’yes’ 

Stay  in  the  program 
yom*  ’  y  ’ 

This  part  of  the  program  produces  the  menu  on  the  screen 

open(4,  file* ’mirror .dat ’ ) 
srita(6,20)  ’[1]  Mirror ’,’[2]  Ho.  of  layers’ 
write(6,25)  borf,HL 
FORMAT (lOX, A,  lOX,  A) 

F0RMAT(14X,A.  IIX,  I) 
print*,  ’  ’ 

VTite(6,30)  ’[3]A1  concentrations: ’, ’high  n  layer’, 

’ low  n  layer ’ ,  ’ spacer  layer ’ 
write(6,36)  xh,xl,xs 
F0RMAT(1X,A,  6X,  A,6X,A,5x,A) 

F0RMAT(30X,f5.2,  12X,  f6.2,llx,f5.2) 
print*,  ’  ’ 

write(6,40)  ’[4]  Design  Havelength(nm) ’ , 

’ [6]  Vary  Layer  Thickness  ’ 


«rite(6,42)  lambdac* 1000, thick 
40  F0RMAT(2X.A,  lOX,  A) 

42  F0RMAT(12X,16.2.  30X,  A) 

Hrite(*,*) 

«rite(6,4S)  * [6]  Oblique  Incidence  (angle)’, 
k  'Polarization  ’ 

writo(6,47)  oblique, '(’ ,anglo, ’)' ,  pol 
48  F0IIMAT(2X,A,  l&X,  A) 

47  F0RMAT(12X,A,,2x,A,f6.1,A,  27X,  A) 

writeC*,*) 
write(*,*) 

HTite(6,S0}  '[7]  Plot  Data  in  mirror.dat’ , ’start  (lun) ’ , 
k  ’end  (nm)’,  ’number  oi  points' 

vrite (6 , 58)  lambdalelOOO , Iambda2el000 ,nla]nbda 
80  F0RMAT(1X,A,  8X,  A,2X,A,2z,A) 

55  F0RMAT(35X,f7.1,  4X,  i7.1,2x,I) 

write(*,e) 
srite(e,e) 

srite(6,60)  ’[8]  Reflectivity  in  the  cavity  or  in  air' 
«rite(6,6S)  reflet 
60  FORMAT (IX, A) 

65  F0RMAT(10X,A) 

Hrite(*,e) 

Brite(*,*) 

srite(6,67)  ’[9]  Maximize  Reflectivity’ 
nrite(6,68)  max 

67  F0RMAT(1X,A) 

68  F0RMAT(10X,A) 
write(*,e) 
write(*',e) 

c 

c  Here  the  interaction  vith  the  user  starts,  allowing  changes 

c  to  parameters  on  the  menu 

c 

print*,  ’Type  in  the  parameter  to  change  (one  parameter 
k  at  a  time,  0  to  proceed’ 
read(*,*) choice 

if ( choice . eq . 1 . and . borf . eq . ’ back ’ ) then 
borf® ’front’ 
nn=l 

elseif (choice. eq. 1. and. borf .eq. ’front’)then 
borf® 'back  ’ 
nn=0 
endif 

if ( choice . eq . 2) then 

srite(*,*)  'number  of  layers  of  the  mirror®?’ 
r*ad(*,*)  nl 
endif 
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if ( choice . eq . 3 ) then 

write(*,*)  'A1  content  of  high-index  layer=?' 
read(*,*)  xh 

HriteC*,*)  'A1  content  of  low-index  layer=?’ 
read(*,*)  xl 

writeC*,*)  ’A1  content  of  spacer’ 
read(*,*)  xa 
endif 

if (choice. eq. 4) then 

write(*,*)  ’design  wavelength  of  mirror  in  nm=?’ 
read(*,*}  lambdac 
lambdac=lambdac* . 001 
endif 

if (choice. eq.S. and. thick. eq. ’no’)then 
thick=’ya8’ 

write(*,*)  ’Thickness  multiplier  for  High  index  layer=?’ 
read(*,*)  factorh 

write(*,*)  ’Thickness  multiplier  for  Low  index  layer=?’ 
read(*,*)  factorl 

elseif (choice . eq.  5 . and . thick. eq. ’yes ’ )then 
thick=’no’ 
factorl=l . 
factorh=l . 
endif 

if (choice . eq . 6 . and . oblique . eq. ’no ’ )then 
oblique* 'yes’ 

write(e,e}  ’Incidence  angle  in  degrees?’ 
read(*,«)  angle 

write(*,e)  'Polarization,  if  TE,  input  1:’ 
read(*,*)  npol 
if (npol.eq.O)then 
pol=’TM' 
else 
pol=’TE’ 
endif 

elseif (choice . eq.  6 . and . oblique . eq. ’yes ’ )then 
oblique* ’no’ 
angle=0 . 
npol=l 
pol=’TE’ 
endif 

print*,  ’angle*  ’, angle 
if ( choice . eq . 7 ) then 

write(*,*)  ’lambdal  (beginning  wavelength;  tuo)  ’ 
read(*,*)  lambdal 
lambdal*lambdal* . 001 

write(*,*)  ’lambdal  (ending  wavelength:  nm}’ 

read(*,*}  lambda2 

Iambda2*lambda2* .001 

write(*,*)  ’nlambda  (number  of  points)’ 

read(*,*)  nlambda 


All 


aadil 

if ( choica . aq . 8 . and . raf let . aq. ' cavity ' )then 
rallct=»air’ 
nn=l 

alsaif (choica.aq.S.and.raflct .aq. ’air’} then 
reflct=’ cavity’ 
il(borf .aq. ’back’)nn=0 
andif 

if (choice. aq. 9. and. max. aq. ’yes’) than 
inax=’no’ 

alsaif (choice. aq. 9. and. max. aq. ’no’)then 
m:uc='yas’ 
andif 

if (choice. na.O) goto  10 
c 

c  initializing  constants  for  this  run 


PI=3.141S92 

angla-anglaaPI/ 180 . 

sina3sin(angla) 

conv=6 . 626a-34a3 . al4/l . 6a-19 

c  datarmina  the  sidth  of  high-indax  (indaxB)  layer  (sidthH)  and  los-index 
c  (indaxL)  layer  (sidthL).  lamdac  is  the  design  frequency  of  the  etalon. 

hvc^conv/lambdac 

CALL  REF(xl,  hvc.  indaxcL) 

CALL  REF(xb,  hvc,  indaxcH) 

top=xh 

nsubs3 . 6 

if (nn. aq. l}nsub=l . 

c  To  provide  max  reflectivity,  ua  need  to  use  the  max  index  changes  at 
c  the  substrata  and  spacer,  or  air  interfaces 
if (max.aq. 'yas')than 

if (mod(nl,2) .aq.O.and.abs(nsub-indaxcl) .It.abs 
A  (nsub-indaxcH))then 

topsxl 
xl=xh 
xh=top 

print*, ’Ent^ing  uith  high  index  to  maximize  reflection’ 
alsaif (mod(nl, 2) .aq. 1 .and.abs(nsub-indexcH) .It.abs 
A  (n8ub-indexcL))than 

top=xl 
xl^xh 
xh»top 

print*, ’Ending  with  low  index  to  maximize  reflection’ 
andif 
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Brite(*,*) 

ii ( top-xa . eq . 0 . ) then 
print^'.  'layer  next  to  cavity  adds  no 
t  reflectivity,  discard?' 
read(*,*),  yorn 
if  (yom.ne.  'n')then 
nl=nl-l 
top=xl 
xl=xh 
xh=top 

print*, 'number  of  layers  is  now:  ',  nl 
print*,  'layer  next  to  cavity  has  x 
k  concentration;  ',  xh 
ariteC*,*) 
ariteC*,*) 
endif 
endif 
endif 


CALL  R£F(xl,  hvc,  indexcL) 

CALL  REF(xh,  hvc,  indexcH) 
if(nn.eq.O)  call  ref(0.,  hvc,  nspace) 
sidthB  =  lambdac/4./indexcH*FACTQRH 
widthL  *  lambdac/4./indexcL*FACT0RL 

if (xl.gt.xh)then 


print* , 

"widthH  =", 

WidthH,  " 

IndexcH  =  ", 

indexcH 

print* , 
else 

"widthL  =", 

WidthL,  " 

IndexcL  =  ", 

indexcL 

print*. 

"widthH  =", 

widthl,  " 

IndexcH  =  ", 

indexcl 

print* , 

"WidthL  =", 

widthh,  " 

IndexcL  =  ", 

indexch 

endif 

writeC*,*) 

write(*,*) 

c  calculate  the  transmissivity  trans  of  the  etalon  as  the  wavelength  is 
c  veuried  from  lambdaKum)  to  lambda2(um). 

XX1=0. 

XX2=0. 

PP=0. 


do  k=  1 ,  nlambda■^  1 

lambda=lambdal* (lambda2-lambdal )/f loat (nlambda) *f loat (K-1 ) 

pi  ^bda=  1  ambda/4 /p  i 

hv=conv/lambda 

CALL  REF(xl,  hv,  indexL) 

CALL  REF(xh,  hv,  indexH) 
call  ref(xs,  hv,  nsub) 
call  ref(xa,  hv,  nspace) 
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if (im.eq.l)nsub=l. 

cindexh=cinplx  (  indexh  ,-0.0) 
cindexl=cinplx(indexl,-0.0) 


do  m  =  2,nl,2 

index (m)  =  cindexL 

C2ill  pha( index (m) ,BidthL,cosd(m) . phase (m)) 

end  do 

do  m  =  l,nl,2 

index (m)  =  cindexH 

call  pha( index (n) ,widthH.cosd(m) ,phase(in)) 

end  do 

c  calculate  the  transmission  trans  at  savelengtn  Icunbda. 

c  refer  to  Computer-aided  Techniques  for  the  Design  of  Multi-layer  Filters 
c  (pp.9-10)  by  H.M.  Liddell(1981,  Adam  Hilger). 

if (numbera.eq. 1)  then 
do  m=l,nl 

if (npol . eq. 1)  index(m)=index(m)*co8d(m) 
if(npol  ne.l)  index(m)=index(m)/cosd(m) 
end  do 
endif 

r  =  ccos(phase(l)) 

s  =  (0,1.)  *  cain(phase(l))  /  index(l) 
t  =  (0,1.)  ♦  csin(phase(l))  ♦  index(l) 
u  =  r 


do  m=2,nl 


rr  =  ccos(phase(m)) 

ss  =  (0,1.)  ♦  csin(phase(m))  /  index (m) 
tt  =  (0,1.)  *  csin(phase(m))  *  index (m) 
uu  =  rr 


end  do 
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b=r+s*nsub 
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c=t+u*nsub 

trans (K) =4 . ♦nspace+nsub/ cabs (ni  pace*b+c ) **2*100 . 
REFL=cabs((nspace*b-c)/(nspace*b+c)) **2*100. 

WRITE(4./0)  ’  AHBDA*1000. .REFL,  trans(k) 

if (LAMBDA* . 00 1 . go . lambdac .  amd . lambda- . 001 . le . lambdac) then 
print*,  "Wavelength:  ",  LAMBDA*1000,  "  refl:  ",  REFL 
endif 

end  do 

70  F0RMAT{1X,3F10.5) 

close(4) 
ariteC* ,*) 
arite(*,*) 
print*, ’continue?' 
read(*,*)yorn 
ilCyom.eq. ’y’)goto  10 
EMD 

c 

c  This  subroutine  calculates  the  phase  change  in  passing  thru 
c  each  layer  lor  given  values  of  the  refractive  index  (  xindx  ),  the 
c  thickness  (  aidth  ),  the  incident  angle, 
c 

subrout in e  pha ( x ind ex, aidth,cosd, phase) 
complex  xindex, phase 
common  sina,xlambda 
sinanssina/reaKxindex) 
as inan=as in ( s inan ) 
cosd=cos(asinan) 

phase  =  2.  *  3.141592  *  xindex  *  aidth  *  cosd  /  xlambda 

return 

end 

c 

c  This  subroutine  calculates  refractive  index  of  AlGaAs  for 
c  given  values  of  A1  composition  (  xconc  )  and  aavelength  (  ev,  in 
c  units  of  eV) . 
c 

SUBROUTIHE  REF(Xconc,  ev,  XI) 

aO  =  6.3  *  19.0  *  xconc 
bO  =  9.4  -  10.2  *  xconc 

EO  =  1.426  +  1.155  *  xconc  +  0.37  ♦  xconc  *  xconc 
EOD  =  1.765  +  1.116  *  xconc  +  0.37  ♦  xconc  *  xconc 
xki  =  ev  /  EO 
xkiso  s  ev  /  EOD 

fx  =  (2.  -  sqrt(l.*xki)  -  8qrt(abs(l .-xki)))  /  xki  /  xki 

fxso  =  (2.  -  sqrtd .*xkiso)  -  sqrt(abs(l-xkiso)))  /  xkiso  /  xkiso 

val  =  EO  /EOD 

epis  =  aO  *(  fx  +  0.6  •  val**1.6  •  fxso)  ♦  bO 
xn  =  sqrtCepis) 
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XX  =  3.65  +  (ev  -  1.6)/. 3  *  .13 

if  ((ev.ge.l.4.or.xx.lt.xn) .and.xconc.eq.O.)  xn 

T«tum 

end 


XX 
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Vertical  Cavity  Surface  Emitting  Lasers  (VCSELs)  are  a  form  of  semiconductor  laser  which  have  their  cavity 
oriented  orthogonally  to  the  plane  of  the  wafer.  The  orientation  necessitates  short  cavities,  highly  reflective 
mirrors  and  a  relatively  high  gain/loss  ratio.  Even  so,  the  resultant  superior  exit  beam  characteristics  and  the 
tight  packing  density  of  the  finished  lasers  provide  strong  motivation  for  pursuing  the  growth  of  these  structures. 
This  thesis  details  the  design  of  an  optically  pumped  InGaAs  multiple  quantum  well  periodic  gain  structure 
VCSEL  with  a  950  nm  lasing  wavelength.  These  growths  were  to  be  a  first  attempt  at  VCSEL  construction,  so 
part  of  this  study  included  verification  of  the  quality  of  the  parts  of  the  finished  design.  These  measurements 
required  the  construction  of  a  laboratory  configuration  to  optically  pump  VCSELs  and  characterize  them  by 
spectral  reflectivity,  output  beam  polarization,  mode,  lasing  wavelength,  and  optimal  pump  wavelength.  Analysis 
of  the  characteristics  for  several  VCSELs  obtained  from  the  University  of  Arizona,  and  the  back  mirror  grown 
locally,  illustrate  the  ability  to  use  measured  data  and  theoretical  spectral  reflectivity  calculations  to  determine 
the  quality  of  the  growths. 


Semiconductor  Lasers,  Gallium  Arsenide  Lasers,  Semiconductor  Devices,  Surface 
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